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ABSTRACT 


A new calculational method for multistage, multicomponent 
separation processes is presented and evaluated in this thesis. A 
modified Newton-Raphson procedure is used to solve all the model 
equations simultaneously for corrections to temperatures, flow rates 
and compositions. 

A detailed description of the method proposed by Tomich is 
given and the new method is compared with Tomich's method and with 
Bubble Point and Sum Rates procedures. These comparisons indicate 
that the new method is a significant improvement on current procedures. 

The method is reliable for both distillation and absorber 
type problems, is particularly effective for reboiled absorber 
calculations and readily handles complex column configurations and 
the direct calculation of composition dependent equilibrium and 


enthalpy data. 
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CHAPTER | 


INTRODUCTION 


Multicomponent, multistage separation processes are commonly 
used in the chemical, natural gas, petroleum and related industries 
and the simulation and the design of these processes are often based 
on the steady state solution of the equations describing an equilibrium 
stage model. The model equations can be numerous and are nonlinear 
so that solution procedures are complex and numerical methods must be 
employed. As a consequence, many calculational methods have been 
proposed. 

The traditional approach to designing a procedure for the 
rigorous solution of the model equations has been to avoid the larger 
problem of attempting to solve all equations simultaneously by 
partitioning the problem into smaller sub-problems. This is 
accomplished by decoupling the effects of certain variables and 
solving selected groups of equations in a particular order for 
compositions, temperatures and flow rates. The equations are grouped 
in such a way that solutions of the subsets are readily obtained after 
assuming certain variables. A convergence technique is used to link the 
subsets so that over-all consistency is obtained. Well known procedures 
such as those proposed by Lewis and Matheson (1), Holland (2), Thiele 
and Geddes (3), Wang and Henke (4) and Sujata (5) are in this category. 
Such methods, based on decoupling, have inherent difficulties in that 
the ability to converge will vary widely from application to 
application and there is no assurance that such methods will converge 


except experience with particular classes of problems. 
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A general algorithm that would be reliable for a broad range 
of feed components or column configuration involved, has been desired 
and it has been recognized that a method based on solving al] equations 
simultaneously would probably provide the ultimate solution procedure. 
However this is a formidable problem bags of the computational 
effort involved. 

The Newton-Raphson technique is an effective method for 
solving complicated sets of non-linear algebraic equations and several 
authors have recently reported on the use of Newton-Raphson procedures 
in iterative methods designed to solve the equilibrium stage model 
equations for several or all variables simultaneously. Examples are 
the methods proposed by Tomich (6), Gentry (7), Naphtali and 
Sandholm (8), (9),Goldstein and Stanfield (10), Billingsley and 
Boynton (11), and Tierney and coworkers (12), (13). These methods 
generally require fewer iterations for solution than do the 
traditional methods and exhibit good convergence characteristics for 
a wide range of problems. However they have the disadvantage that 
the large number of partial derivatives that must be evaluated and 
the methods used for solving the linearized sets of equations 
generally result in large computer storage requirements and long 
computational time. The methods may become particularly unwieldy 
when composition dependent enthalpy and equilibrium data are employed. 
In addition several of these methods still decouple the effect of 
composition in that the component material balance equations are 
solved separately. Thus, there is still plenty of scope for devising 
more efficient methods for solving the equilibrium stage model equations. 
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the author developed the bsic theory for a new calculational procedure 
that showed promise of being an improvement on the procedures presently 
available. 

The objectives of this investigation are: 

(i) to elaborate the basic theory for practical use. 

(ii) to prepare a computer program based on the new method 
and evaluate its performance for distillation and 
absorption type problems. 

(iii) to prepare a computer program based on the algorithm 
suggested by Tomich (6), which appears to be one of 
the most promising of the recent methods, and to 
compare its performance with the new algorithm. 
The new method makes use of a linearization technique which 
is actually a multivariate Newton-Raphson procedure. The method 
solves all the linearized equilibrium stage model equations 
simultaneously. Convenient matrix structures are formulated by the 
use of prudent assumptions when selecting the linearized equations 
and these matrices are solved using a procedure for inverting tri- 
diagonal matrices that had previously been developed by the author. 
The method is characterized by the simultaneous solution of the 
linearized equations and by the use of new matrix manipulation pro- 
cedures. In addition the new algorithm takes into account the com- 
position dependency of enthalpy and equilibrium data. 
In the computer programs the Chao-Seader correlation was used 
exclusively. It is representative of generalized correlation procedures 


and it has been widely accepted for use in the hydrocarbon processing 


industries. 
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CHAPTER 2 


LITERATURE REVIEW 


2.1 Equilibrium Stage Model 

A generalized equilibrium stage model is shown as Figure 1. 
This model may be used to represent a wide range of separation devices 
including complex column configurations where multiple feeds, side 
streams and/or side heaters and coolers may be employed. Each stage 
in the model is assumed to be an equilibrium stage, that is, the streams 
leaving a stage are assumed to be in thermodynamic equilibrium. 

The steady state model equations for any stage j and any 
component i may be represented as follows. 


(1) Equilibrium Relations: 
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(5) Over-all Material Balance: 
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The component material balance equations are often combined 
with the equilibrium relations as, 
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There are N x (2 x NCP+3) independent equations and 
N x (2 x NCP+3) unknowns provided that the following quantities are 
specified, 
(i) the number of equilibrium stages 
(ii) the pressure at each stage 
(iii) the rates, compositions and conditions of all feed 
streams 
(iv) the rates of all side streams 
(v) the amount of heat added or removed from each stage 
including any condenser or reboiler. 
The solution of the model equations then requires finding values for 
the temperatures, eto the liquid and vapor flow rates, Eas and 
V.'s, and the phase compositions, x .'S and hes at each 


J J 
Stage such that We Og Ori E and = are all equal to zero. 


lJ 
The specifications listed above are those normally made for 
an absorber type problem. It is permissible and often desirable to 
interchange some of the specified variables with the unknown variables. 


For distillation problems it is usual to specify the amount of top 


product and the reflux ratio and to treat the reboi ler and condenser 
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loads as unknowns. Similarly, the top vapor flow rate, V is 


12 
usually specified for reboiled absorber problems and the reboi ler 


load is treated as an unknown. 


2.2 Solution Methods 

In general, the equilibrium stage model equations form a 
quite complicated and nonlinear system so that numerical techniques are 
required in order to obtain solutions. One approach, and perhaps the 
most mathematically desirable approach, to solving the equations is to 
linearize all of the model equations and to solve them simultaneously. 
Successive solutions would be obtained until a satisfactory convergence 
is achieved. This is referred to as a Newton-Raphson procedure. How- 
ever, the general Newton-Raphson procedure requires the evaluation of 
a large number of partial derivatives and a prohibitively large amount 
of computer storage and calculation time. Consequently considerable 
effort has been directed toward finding more tractable methods. 

Friday and Smith (14) do an excellent job of analyzing and 
categorizing traditional calculation methods that are based on 
decoupling the model equations and clearly outline six decisions that 
must be made to formulate a solution procedure. Published methods 
differ because of the different approaches taken in making these 
decisions. 

The first decision is to choose one of two formulation 
procedures which concern the grouping of the model equations, namely 
grouping by stage or by type. Greenstadt et all (15) and Bergamini (16) 
group the equations by stage. However most methods group the equations 


by type. Generally grouping by type is preferable because grouping 
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by stage essentially makes the method susceptible to build up of 
truncation error. 

The second decision involves the order of satisfaction of the 
equations. The order (1-2)-4-3 is most common. Equations (2-1) and 
. (2-2) are always combined and satisfied simultaneously. Also 


equation (2-5) is used to eliminate one set of the phase rates, uo 


The third decision concerns the selection of the appropriate 
type of equation to provide a given variable. Any set of compositions, 
which are obtained by simultaneous solution of the material balance 
equations and equilibrium relationships, can be used to estimate the 
new flow rates through the summation of the unnormalized compositions, 
or they can be used to supply the new temperatures. The former 
procedure is called the Sum Rates method and it has been applied to 
absorber problems. The latter is called the Bubble Point method. 

Most methods use the latter procedure to obtain the new temperatures 
and then energy balances to generate the new flow rates. 

The fourth decision is the selection of a method to solve 
the material balance equations for the new stage compositions. One 
is the stage-to-stage method which usually calculates from both ends 
of a column toward the middle. The second makes use of a matrix 
which describes all material balance equations simultaneously for each 
component. Stage-to-stage methods are often numerically unstable because 
of the build up of truncation error. This is especially serious for a 
system which contains components whose relative volatilities differ 
widely in magnitude. 


The Lewis-Matheson (1) and the Thiele-Geddes (3) method are 
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representative of stage-to-stage procedures. In the former method the 
feed, reflux and the separation of two key components are specified 
and the number of theoretical stages required for a desired separation 
is calculated. The feed plate is determined by comparing the 
composition obtained by calculating from both the top and the bottom 
of the column. The mismatch at the feed stage is used to correct the 
assumed end compositions. The Lewis-Matheson method is advantageous 
for design work since it computes directly the required number of 
stages for a specified separation and also determines the optimum 
point for feed introduction. However it is difficult to obtain the 
solution for a complex column where multiple feeds or side streams 
are involved. In the Thiele-Geddes method the number of the 
theoretical stages in each section of the column, feed flow rate, 
reflux rate and the quantity of the distillate must be specified. 
The two methods differ in that the Thiele-Geddes method initially 
calculates the ratio of the stage composition to end composition 
rather than stage composition. Then using the material balance 
equation, the ratio of composition at the top plate to that at the 
bottom for each component is calculated. The flow rate of top product, 
which can be calculated from the above ratios, is compared with the 
specified one. If matching is not obtained, the same procedure is 
performed successively. Holland and coworkers (17), (2) have developed 
a convergence technique called ''the Theta method'' which they have used 
with the Thiele-Geddes method to improve the performance of the 
Thiele-Geddes method and to extend its range of application. 

Matrix methods have been used widely. They are not 


susceptible to truncation error and they work equally well for any 
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number of feeds and side streams, and handle nondistributed components 
in the same manner as distributed ones. In these methods specifications 
are the same as those in the Thiele-Geddes method. The first rigorous 
computational procedure for the simultaneous solution of the material 
balance equations was proposed by Amundson and Pontinen (18). Later 

the tridiagonal matrix method incorporating Thomas' algorithm was 
proposed by Wang and Henke (4). 

The fifth decision concerns the stage temperatures. They 
are usually determined by bubble point or dew point calculations. In 
these calculations, the Regula-Falsi method or the Newton-Raphson 
method have been often used. Wang and Henke (4) recommend the Muller 
method for bubble point or dew point calcualtions. However it is quite 
time-consuming to do those iterative calcualtions. To avoid bubble 
or dew point calcualtions, Holland and coworkers (19) developed the 
KE method which uses relative volatilities to a reference component. 
Newman (20) has proposed the use of the Newton-Raphson method for 
solving the stage temperatures so as to satisfy all material balance 
and equilibrium relationships. In this case Thomas' tridiagonal 
matrix algorithm can be used. A similar concept has been exclusively 
used in the Sum Rates procedures (5), (21), (22) to calculate 
temperatures so as to satisfy the energy balance equations. 

The sixth decision is the method of obtaining the new flow 
rates. In the Sum Rates method the new flow rates are immediately 
determined from the calculated compositions. In the Bubble Point 
method the flow rates are determined so as to satisfy the energy balance 
equations. 


The above six decisions deal with the calculation at each 
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iteration. Solutions are attained by successive approximation until 
a certain convergence criterion is satisfied. The following sequence 
is usually employed. 

(1) Assume a set of temperatures and flow rates. 

(2) Calculate phase compositions. 

(3) Correct flow rates and temperatures. 

(4) Return to step (2). 

However, there is no assurance that these decoupled me thods 
will converge. For example, computational procedures which are very 
effective for a close boiling point distillation may have difficulty 
converging for a wide boiling point feed and procedures that work 
well for absorber and extration problems often fail on ordinary close 
boiling point distillation problems. 

For a mixture which consists of extremely volatile and 
extremely nonvolatile components, phase flow rates may be completely 
dominated by the relative volatilities and may not be affected 
relatively by temperatures. If the temperatures are determined by the 
Bubble Point method in this case, a slight change in compositions will 
lead to much temperature fluctuation. Then the ordinary Bubble Point 
method is not suitable for solving this type of problem. This is 
often the case for absorber problems. 

Sum Rates procedures have been proposed for such problems. 
McNeese (21), Sujata (5), Burningham and Otto (22), and Friday and 
Smith (14) presented methods which successfully solve absorber type 
problems. These methods are based on the Sum Rates procedure as 
described in the third decision. However the methods usually cannot 


solve problems dealing with close boiling point mixtures. In addition 
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the Sum Rate methods do not appear to work well for reboiled absorber 
problems (22), (23). 

Holland and coworkers (2),(17) have utilized exclusively the 
Thiele-Geddes method in conjunction with the theta convergence 
procedure to solve absorber and reboiled absorber problems. However, 
in order to obtain convergent solutions of these problems, they found 
it is necessary to employ special techniques to handle the energy 
balances. These are the Q-method and the constant composition method. 
The Q-method introduces side heaters or coolers on each stage so as 
to maintain specified flow rates. However, it does not necessarily 
give a practical solution. The constant composition method of 
calculating stream enthalpy is introduced to avoid round-off errors 
associated with conventional enthalpy balances. This method cannot 
be applied correctly for systems utilizing composition dependent 
enthalpy data. However, their method has been one of the successful 
procedures for solving reboiled absorber problems. Another possibility 
to solve reboiled absorber problems may be to make use of relaxation 
methods (24), (25), (26). Prowse and Johnson (27) have successful ly 
employed a relaxation procedure to solve a reboiled absorber problem. 
The characteristics of these methods will be described later. 

Friday and Smith in their analysis pointed out that an 
improvement in reliability could probably be obtained by solving the 
C-matrix equations (2-6) for phase compositions and then solving the 
E, and $,. equations, (2-3) and (2-4), simultaneously for corrections 
to temperatures and flow rates. Algorithms based on this approach 
have been presented by Tierney and Bruno (12), Tierney and Yanoski C3 )e 


and Tomich (6). Their methods make use of the Newton-Raphson procedure. 
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Namely, functions which show material imbalance or energy imbalance 

are approximated by expanding as functions of flow rates and temperatures 
using a Taylor series expansion truncating after the first order 
derivatives. These functions are then set equal to zero and solved 

to obtain corrections to temperatures and flow rates simultaneously. 

In this case the Jacobian matrix must be evaluated numerically rather 
than analytically. 

Tomich (6) applied Broyden's procedure (28), which is actually 
an improved Newton-Raphson method, to obtain corrections for temperatures 
and flow rates simultaeously. Broyden's procedure updates the inverse 
of the approximate Jacobian matrix using residuals calculated at the 
preceding iteration level instead of evaluating a large number of 
partial derivatives numerically and inverting the rigorous Jacobian 
matrix at each iteration. Tomich's method reduces the amount of 
computation required significantly in that only one matrix inversion 
is required per problem solution. 

Billingsley and Boynton (11) suggest another method to 
reduce the amount of computational effort in the Newton-Raphson 
procedure. They assume temporarily that the material imbalance 
equations are functions of only temperatures and then solve for 
temperatures by the Newton-Raphson method as proposed by Newman (20). 
These temperatures are then substituted into the original equations 
which are generated from the material imbalances and energy imbalances 
by the Newton-Raphson procedure. However, this decoupling assumption 
seems to be often serious, especially for a wide boiling point range 


mixture. 


Friday and Smith also suggested that the ultimate method 
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could be to linearize all the model equations first and then solve 
simultaneously for all variables. Even though this concept has been 
recognized as desirable, it has been difficult to find a practical 
algorithm because of the inherent computational effort. Naphtali and 
Sandholm (8), (9) have proposed a method based on this concept of 
linearization and solving all variables simultaneously. In their 
method all equations are grouped according to stage rather than type. 
The method may take composition dependent data into account. A 
similar method has been proposed by Gentry (7). However these new 
methods are still accompanied by several difficulties. Namely, a 
large amount of computer storage and a tremendous amount of calcula- 
tional effort for the evaluation of partial derivatives are required. 
These problems have restricted the exact application of the methods. 
For example Gentry (7) solved only problems with composition indepen- 
dent data. To avoid such restriction of computer storage Goldstein 
and Stanfield (10) present a new algorithm. In their method enthalpy 
and equilibrium relationships are assumed to be functions of temper- 
atures only and all equations are grouped by type and then linearized. 
When a problem involves a relatively large number of theoretical 
stages, the method makes use of a special technique to cope with the 
limitation of computer storage. For such a problem, a sectioning 
procedure, which calculates for a group of stages assuming linearity 
of variables in the group, is applied. 

Other representative solution procedures are the so called 
"relaxation methods!! and the methods based on an analytical approach. 

"Relaxation methods'', in principle, calculate the transient 


behavior of the stage variables until steady state operation is 
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approached. Flow rates, stage compositions and temperatures are 
first assumed and then the variables corresponding to ene ae are 
calculated so as to compensate imbalances in enthalpy and material 
balance at each stage. These procedures are continued until all 
equations are satisfied. Rose, Sweeny and Schrodt (24) have 
successfully developed this approach. Applications and improvements 
have been proposed by Hanson, Duffin and Somerville (29), Ball (25), 
Ishikawa and Hirata (26), and Prowse and Johnson (27). Relaxation 
methods are highly stable numerically. Their high stability can be 
helpful when the system contains a wide boiling point range mixture 
or a highly non-ideal mixture. However, the rate of convergence of 
relaxation methods is quite slow, especially as the solution is 
approached. 

Acrivos and Amundson (30) present the analytical solution 
of an ideal distillation problem for the case when constant mclar 
overflow and constant relative volatility can be assumed. Acrivos 
and Amundson (31) utilized perturbation techniques to extend the method 
to non-ideal systems. This idea has been further extended by combining 
the concept of the theta convergence procedure by Yamada and Sugie (32). 
However the methods for a non-ideal systemare quite complicated and 
it does not seem to be practical. Furthermore it is not easy to take 


account of energy calculations to generate new flow rates 
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CHAPTER 3 


THE TOMICH METHOD 


The procedure proposed by Tomich (6) for solving the 
equilibrium stage model appears to be one of the most promising of 
the recent methods. Thus it was selected as an appropriate method 
to use as a basis of comparison when evaluating the performance of 
the calculational procedure proposed by the author. In Tomich's 
method the effect of composition is decoupled in that the component 
material balances are solved separately. Corrections to temperatures 
and flow rates are then determined simultaneously using a multivariate 
Newton scheme. Broyden's procedure is effectively applied to reduce 
the computational effort required for the evaluation of partial 
derivatives. Details of this method follow. 

The component material balance equations combined with the 
equilibrium relationships, that is equations (2-6) with C. ein OF 


’ 


may be expressed in a tridiagonal matrix form. 
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For given liquid and vapor nee and equilibrium ratios 
these equations are a linear set which may readily be solved by using 
the Thomas algorithm for a tridiagonal matrix as described by Wang 
and Henke (4). 

Imbalances for energy and the summations of compositions 


at each stage are given by, 
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In the Tomich method the temperatures, ie the flow rates, 
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Vos eeule. Se) andsthe compositions, ox. -'s y. .'S, are calculated 
J J I,J I,J 


by using the Newton-Raphson method so as to satisfy the following 


conditions: 


E,=0, (3-4) 


S; =0. (3-5) 


A multivariable nonlinear function may be approximated in a 
linearized form by expanding the function as a Taylor series truncated 


after the first order derivatives. That is, 


k+] k 
f (xX) sXpo+++ 9%) ~ f (X) 2X5 oe+5 9X) + \ a Ax (3-6) 


where k is the iteration number and the x Ss are variables. 


=f and § are assumed to be functions of the temperatures 


and the vapor flow rates only, and are approximated as follows: 
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Since the * and * are errors or residuals in the 
summation and heat balance equations and are equal to zero when a 


solution is reached, these relations are set equal to zero at the 


‘ k+] 
new iteration level. The equation set (3-7) and (3-8) with EB = 0 


+ : , d 
and 5 = 0 can be written in a matrix form, 
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(3-9-a) 


The strict application of equation (3-9-a) is to an absorber 


the flow rates, 


Ve 
J 


's are unknowns. 


For 


normal reboiled absorber and distillation type problems slight 


modifications are required. 


are as follows, 
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For a reboiled absorber where 
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Matrix expressions for these problems 
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V is specified. 


(3-9-b) 
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For a distillation problem where V, and V, are specified. 
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From equation (3-9-a-b-c) the corrections, ee and 


a at are calculated. Then ur s and se are corrected as 


follows, 
7 OT EAY (3-10) 
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where t is a weighting factor. The determination of t will be 
shown later. 

In equations (3-9-a-b-c) E and Sj are calculated from 
the definitions (3-2) and (3-3). The partial derivatives in the 
Jacobian matrix must be evaluated by a finite difference method. 

The eer ee sets (3-1) and (3-9-a-b-c) are then solved and the same 
calculations are performed iteratively. The basic procedure may 


be as follows, 
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(1) Some initial set of Vits and Tits is assumed. 

(2) The L.'s are computed from equation (2-5) with M =a Oe 

(3) The equation set (3-1) is solved by using the tridiagonal 
matrix algorithm to give the A esas 

(4) The y. .'s are calculated from equation (2-1). 


les 

(5) The E's and Si's, and then partial derivatives are 
evaluated. 

(6) Equation (3-9-a-b-c) is solved for AV."s and AT; 's. 

(7)5, Ihe Vi"s and Wis are calculated from equations (3-10) 
and (3-11). 

(8) By using this new set of Vits and Tits, steps (2) 
through (7) are repeated until an appropriate convergence 
criterion is satisfied. 

However, the procedure mentioned requires the evaluation 
of the Jacobian matrix numerically at each iteration. In this case 
it is quite time-consuming to compute all partial derivatives in 
the Jacobian. Furthermore a large computer effort is required to 
obtain the inversion of the Jacobian matrix. To cope with these 
difficulties Tomich applied the Broyden procedure (28). The 
Broyden procedure is a modified version of the Newton method and 
is designed to reduce the number of function evaluations required and 
to update the inverse of the approximate Jacobian matrix. An outline 
of this procedure follows. 


Consider the set of n nonlinear equations, 
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These may be written more concisely as, 
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f(x) = 0. (3-13) 


sKee. : 
fe x is the cee approximation of the solution of 
equation (3-13), then the next variables are defined as follows 


according to Newton's method: 


eee ene (3-14) 


where a‘ is the Jacobian matrix. 


Let aK be defined as 


Pa -B c 9 (3-15) 


wk : ; : 2 
where B is some approximation to the Jacobian matrix A. 


Then a simple modification to Newton's algorithm gives 


by 
giten= 35 48r5c8 (3-16) 
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= (3=17) 
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B is chosen so that the change in f predicted by 8B 


in a direction orthogonal to ae is the same as would be predicted by 


BK. Then the next relation is obtained by using Householder's formula. 
T 
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where t is a weighting factor and p is a transpose vector of p. 


Be is actually an approximate inverse form of the Jacobian 


matrix ma Then the new x is calculated as, 
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= x + tebe (3-20) 


The determination of a weighting factor t, which moves the 


solution toward convergence, is described by Broyden as follows. 


The first value of a weighting factor 2) at each 


iteration is chosen to be unity since this is the value arising 
naturally in Newton's method. The second value, if required, is 


given by the semi-empirical relation, 


(2). (+60)? - 1 


xe (3-21) 


where 


ae ah (3-22) 


and where $(t) is the square of the Euclidean norm of f(t). 
lf further improvement to ¢ is required, it is approximated 
by a quadratic function whose ordinates are specified at t = 0,1 


(2). (3) 


and If this function is convex, then is chosen to be 


the value of t that minimize $4, and a new quadratic is formed 


Q) 


whose ordinates are specified at and two of the previous trial 
values of t. If there is no value of t, (0 < t <1), such that 
o(t) is smaller than 9(0), then t will be searched for between 
DeeoliGame— 1. 

Functions, ae in the general equation (3-12) represent 
the residuals in the model equations, namely, Aye and Spee The 
Vector. x, in the general equation corresponds to the temperatures, 
RSE and the vapor flow rates, eu in the model equations. 


Consequently the approximated inverse of the Jacobian matrix, which 


is calcualted from equation (3-19), is expressed in terms of yee 


Hoes 36 w, 10934) gotadylaw & 10: s0tew Jed RT | 
Gitlehew ute ‘ars et gina ean e\qd Ina.ad oF noveta af notreved! 
et tenth tips 7i ,adlev bnedse sat Jbondem <!notwel al yitensen 
folselen best 1 1qma- lange ett ye nevi - 


(it+e) La : 
(se) . ‘" 


4 Io mon assbi {aug az io sjevpe ody zr (306 sreite bee 


bededtxoigoe 2) 2) .oorlupes eb @ 0% Srsmevcnan) sot ail FI - 

PO = 4 de baliigage tye esisnibro seordw noisoent alievhaup & yd | 

$d.62 needs 2i El, ney .xsvauo ei cobtand? 2ida ¥4 mM, bre : 
hertiot af ofserbeup won s bre .f aximioin deni 1 eo wiley ata’ 
feted eumivarg ed Yo ow? tne hy ja bellisege ais 2stecliw egote 

we, Cos a 2 4 > 0) 3 Tolnutoy on ef sraty FI 4% eoulev | 


pentane at ht +oe an asi? — fon, 


24 


Seliss Rds and Nhs at present and previous iteration levels. 

— Now for succeeding iterations, Broyden's procedure is 
effectively used for updating the inverse of the Jacobian matrix 
using information at preceding iteration levels. Thus there is no 
need for more than one evaluation of partial derivatives and matrix 
inversion per problem solution and the computation time may be 
greatly reduced. 


The improved calculation procedure may be summarized as 


bl) edSomel initialhvseteof Vils and Tits is assumed. 

(2)4 The Li's are computed from equation (2-5) with oF = 0. 

(3) The equation set (3-1) is solved by using the tridiagonal 
matrix algorithm to give the ieee 

(4) The Vij s are calculated from equation (2-1). 

(5) The E i's and S's are calculated by substituting the 
above information into equations (3-4) and (3-5). 

(6) The approximate inverse of the Jacobian matrix is calculated 
based on the Broyden procedure and the new variables are 
calculated from equation (3-20). 

(7) By using this new set of V's and Tits, steps (2) 
through (6) are repeated until an appropriate convergence 
criterion is satisfied. 

However, at the first iteration only, all the partial 
derivatives must be evaluated and the inversion of the Jacobian matrix 
must be obtained. The procedure for doing this is the same as the 


general Newton-Raphson method. The evaluations of the partial 


derivatives are performed numerically, because the aye in 


ions 


a. 


on zt: svarnl eunT siadat ae 2 


even en ait os E791 


ms 


inten fina #ev boa) 9b fe (45g 40 he ic sdutevs sno rei 9 

5d uit auls Ol s67ugneD oil? bnk, didlo eain ae 

auae 

ge Westremmue od yan ewwhsoorqreolteliglss bavosget ST 

.bsmedes 21 ea! bins i qo, Jee Isto ini’ smoe (#) 

0 = .M Ri lw (2+S) noiteups moy) basugnoe sin <4) 1 AF (t}-& 

fenoos tb fay ods unlaw va Kavice ei (lE) S52: policops stT (a) _— 
ste ans _ayip°o2. meittopis xen >) 

(ts) ited mort haselariss 616 hy onl 

otf onlsuyioaduz, ve batelupts> sre at bne a I ant 

¢(2-8) bos UN-E) ena ljeups ain) nolieniolhi aveds 1 

besalustes €) xivicm refdosel ay Fe e2iaynl sJamixortdqe SAT oO 

S16 astdelsey wan siz baa. ayuhssrig nshyor8 sil oo. beeed 

‘ (os- ~E) ‘no Lieype ment bear uotieg 

(S}-2qgsa na bne 2' W. te ae win ely poten a eo 

saileg eves ZiGinqerngs n6. ne batesasy ne (6). ipwonty 


~~ - 


_ sbaitedoes “i adil 


— ia sre nea) a ext aan wid Fe 


‘zh ahs 


oe 


25 


equation (3-2) and the S's in equation (3-3), which show the 
residuals in the energy balance equations and the component material 
balance equations, cannot be expressed analytically. If K-values 
and enthalpies which are dependent on compositions are employed 
initial phase compositions must be assumed at step (1). The initial 
assumptions are often provided by substituting ideal K-values into 
equation (2-6) and solving the C-matrix equations with Cia 0 by 
the Thomas algorithm. 

In the computer program developed calcualtions are 


continuied until a convergence criterion such as the following is 


satisfied, 


{eR 0°® (3-23) 


where [CRIT] is defined as follows according to Tomich (6), 


N 1 
[CRIT] = J (s5 + Ey) (3-24) 
JF] 


and where 
pa ee eter ts (3-25) 
VaR eH ETLOF a be HheEpeaeie haters 
In equation (3-25), E is normalized by dividing by the 
total heat input to any stage j so that it is of the order of 
magnitude 1. $, the square of the Euclidean norm of f, in 
equation (3-22) corresponds to [CRIT] in the program. Then a 


weighting factor, t, is searched for to satisfy the following 


relation according to Broyden (28) as stated before. 


tccriti<*! «< cerity* . (3-26) 
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CHAPTER 4 


EQUILIBRIUM AND ENTHALPY DATA 


In the range of low or moderate pressures remote from the 
critical conditions, equilibrium K-ratios can be closely approximated 
by using the vapor pressure of each component and the system pressure. 
For a non-ideal system K-ratios are often modified by liquid phase 
activity coefficients. The Wilson correlation (33) is typical of 
such a modification. 

Phase equilibria at high pressures such as encountered in 
petroleum or natural gas processing industries are more complicated. 
The convergence pressure concept has been used commonly to correlate 
K-ratios in a multicomponent hydrocarbon system at high pressure. 

The N.G.P.S.A. charts (34) are the most common source for the 
prediction of K-ratios. The B.W.R. equation (35), (36), which is one 
of the first analytical methods to correlate composition dependency, 
has been often used to predict equilibrium K-ratios. A general cor- 
relation for calculating equilibria in hydrocarbon mixtures was 

first proposed by Chao and Seader (37) and has been widely accepted 
in the hydrocarbon processing and related industries. Attempts to 
elaborate a general correlation are continuing. The Prausnitz and 
Chueh method (38), (39), and the Lee and Edmister method (40) are 
typical of those. 

The enthalpy of a pure component is often expressed as a 
function of the temperature and pressure of the system and the 
enthalpy of a mixture can be evaluated from these pure component 
enthalpies. Yen and Alexander (41) have proposed general enthalpy 
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correlations based on the theory of corresponding state. Edmister et 
al (42) have developed procedures for enthalpy calculations based on 


thermodynamic relationships and the Chao-Seader correlation. 


4.1 Chao-Seader K-ratio Correlation 

The Chao-Seader correlation has been widely accepted in 
hydrocarbon processing industries. In this thesis the correlation 
is exclusively employed to generate equilibrium and enthalpy data. 

Physical properties used with the Chao-Seader correlation 
to calculate equilibrium and enthalpy data are shown in Appendix 1. 

This correlation utilizes vapor fugacity coefficients, 
liquid fagacity coefficients and liquid activity coefficients. The 
vapor fugacity coefficients are calculated by the Redlich-Kwong 
equation of state (43). The liquid fugacity coefficients are 
evaluated by the curve fitted relationship and the acentric factor 
concept of Pitzer. The solubility parameter and the theory of 


regular solutions are applied for calculating the liquid activity 


coefficients. In the correlation there are no parameters which depend 


on interactions among the molecules in a multicomponent system. 
The equilibrium K-ratios may be defined as follows by 
setting the fugacity of any component i in the vapor phase equal 


to that in the liquid. 


eG Nt (4-1) 


Chao and Seader correlated the non-ideality of each 
component in both the liquid and vapor phase to the properties of 


pure components. The vapor phase fugacity coefficient, Vis is 
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calculated based on the Redlich-Kwong equation of state. The liquid 
phase activity coefficient, Yio is calculated from the Scatchard- 
Hildebrand equation for regular solutions. The liquid fugacity 

coefficient, v; = ede 


pressure, reduced temperature and the acentric factor. 


/P, is curvefitted in terms of reduced 


4.1.1 Evaluation of Vapor Fugacity Coefficient v 
In the Chao-Seader correlation the basic equation for 
calculating the vapor phase properties is the Redlich-Kwong equation 


of state (43), namely, 
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where Z_ is the compressibility factor, and A and 8B are constants 


expressed in terms of the critical properties of a component as, 
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The fugacity coefficient vi; is commonly expressed as, 
chotal ter oe }2P PhRT he > 
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Oi 4 Vas 
| J 


; y" Jypiat Feed, yaigenud, ; wast tn sobsmul en ol wlat : 
Tel WHINE feed ait nol telemnes Astestepend otha! © 7M 
Aetieups wirow-nolMbAs sit) el est inagotg seatg rousy wiht potiefustes 7 ; 


elses (ON) amore Ro 
£ 
ct | 
-t a” eee ee, 
alia | ch = Ma - 
tas") = = 7A = ara : 
of 


dinetenos oe 8 oboe A brs Soryet yil fictions yqmen itt ef 3 edw 
fe JNQHagOs = 1p eaitasqate (ealsio si? Yo ema al bere we. 


tene-#) 


fd>e-") 


> oe 


29 


This equation is combined with the Redlich-Kwong equation 
and then the fugacity coefficient of any component i in the vapor 


phase is obtained, 


ai a? a a BP 
Ln, = Zals 7: 2n(Z-BP) a Re ee 3 a) en a aa, (4-5) 
4.1.2 Evaluation of Activity Coefficient Y; 
The liquid phase activity coefficient Y; is calculated from 
the Hildebrand and Scatchard equation, 
[iA ie é) 
any spies ig esac (4-6) 
AN. 
65 GH (4-7-a) 
Vv. 
’ 6.xv 
c= — (4-7-b) 
) xv 
where ve is the molal volume of the liquid phase of component i, 
AH is the heat of vaporization and oF is the solubility parameter. 
The equation for calculating Y; is based on the assumption 


of the regular solution for which the excess entropy of mixing is zero. 
The Gibbs excess free energy consists of enthalpy and entropy of 
mixing. Therefore the non-ideality in a regular solution is entirely 
due to its heat of mixing only. Generally, the assumption of the 
regular solution is not valid. However, systems of non-polar fluids 
such as hydrocarbon mixtures behave close to the regular solution 
theory over a relatively large range of pressure, temperature and 
composition. The characteristic of the method is that the evaluation 


of the activity coefficient Yj is performed from the assumption that 
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Yi; is a function only of the properties of the individual pure 


components. 


4.1.3 Evaluation of Liquid Fugacity Coefficient v, 

The theory of corresponding states has been extended by 
Pitzer introducing the third parameter W; which is called an acentric 
factor. Any property of a fluid is then assumed to be given as 
function of reduced pressure, reduced temperature and the acentric 
factor. 

Chao and Seader applied this concept to correlate the 


liquid fugacity coefficient, Vir as, 


log Sed log fe + w, log ie (4-8) 


Ao is v. of a simple fluid state and ee accounts for 


where 


the departure of v; from that of a simple fluid state. 

However, when the pure liquid of any component i of the 
mixture becomes hypothetical at a given system condition, v; cannot 
be evaluated. In order to overcome this difficulty they extended the 
liquid fugacity coefficient correlation into the hypothetical region 
using experimental equilibrium data. Then they curvefitted the 
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expression Vi; and with approximating functions of reduced 


properties, 
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(0) 2 3 
log Vi = A + Ao/Tp + A3T. - A,Tp + AcTp 
2 
- (A, + ADT, + AgTR)P. 
2 
sal Oakey td Wr Sop deal ete, (4-9) 
Gh) 3 
log oF * Ay, + Aral. aw Aj3/Tp ate Aryl se Ay, (Pp-0.6) (4-10) 


where all the constants A, whichare common for all components, were 
obtained for frequently occurring hydrocarbons. 

The liquid fugacity coefficient for any component i is 
then calculated only from the reduced pressure, the reduced temperature 


and the acentric factor using equations (4-8), (4-9) and (4-10). 


4.2 Enthalpy Calculation From the Chao-Seader Correlation 

The general correlation of enthalpy data based on the 
thermodynamic relationships and the Chao-Seader correlation has 
been developed by Edmister, Persy and Erbar (42), and Erbar (44). 


The ideal partial molar heat of vaporization is given by, 


——=—-) (4-11) 


The enthalpy of mixing due to the non-ideal solution, 


assuming the excess entropy of mixing is zero is obtained from, 


ie 9 a) SOR Y; 
h, -h, =- —— : 4-12 
A 15 RT” ( oT phi ala) 
Then the partial molar enthalpy of component i in the 


liquid is obtained from equations (4-11) and (4-12). 
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where the partial derivatives as obtained from the Chao-Seader 


correlation are, 


es 2 
) aS (6,-8 ) (4-14) 
oT year. 2 ion 
R he RT AT. 
d2n Vv; A, , 
——— = 2.303 {- —+ 
( aT, dp. 303 o ~ Azt2A,TR+3A,T, + (AD+2AgTp) PE 
R 
A 
2 Le) 2 
eo oe ine (4-15) 
T 
R 
Then molar enthalpy of the liquid mixture is, 
NCP . 
rersaq Duel Dinh wpe. (4-16) 
i=] 
The deviation of the molar enthalpy of the gas phase is 
given by, 
H HS fe tv-T(24) dp. (4-17) 
- oT p 


The molar enthalpy of the gas phase mixture is obtained from 
equation (4-17) and the Redlich-Kwong equation of state. 


CP 2 
H= ) y HO - 2 on) + SP) + 1-zdRT (4-18) 


The molar enthalpy of component i as an ideal gase, ee 


is evaluated from the cubic equation in temperature °F. 


He = (AH) . + (BH) .T + (cH) .T° + (DH) .7° : (4-19) 


4.3 Hypothetical Components 


Complex hydrocarbon mixtures such as an absorber oil consist 


of components of wide boiling point range. It is highly desirable to 
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include the Nii ee 48 handle complex hydrocarbon fractions as well 
as pure components. Erbar (44) and Cavett (45) have developed 
procedures to predict the properties of complex mixtures to provide 
data for the equilibrium and enthalpy calculations. The minimum 
information required for the calculations is the API gravity, the 
molecular weight and either the mean average boiling point or the 
molal average boiling point and the cubic average boiling point of 
the fraction. From this minimum information the equilibrium and 


enthalpy calculations are performed as, 


(1) If the mean boiling point is not specified it is calculated 


from the molal average boiling point and the cubic average 


boiling point by, 
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pp = —ca__mba (4-20) 


(2) “Ghe specific arav ity of ithe fract tongat — 60°4P, is 


calculated by means of, 


specific gravity = oo (4-21) 
: AP | 


and the density is adjusted by 25° C by, 


d,- = 0.98907 x specific.gravity at 60° F. (4-22) 


25 


(3) The liquid molal volume at 25° C_ is calculated by, 


* — molecular weight (4-23) 


V = 
dog 


(4) The pseudo-critical temperature of the hypothetical component 


is calculated by, 
2 nag 2s? 
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where T is the molal average boiling point, if available, 
or the mean average boiling point in °R, and A is API 
gravity. 

The pseudo-critical pressure (psia) of the hypothetical 


component is calculated from, 


~ 2 3 2 
log Pe = bytb  T+boT el +b, AT+b (AT 


2 219 
oA T+b A hea (4-25) 


+b 
If the critical pressure has been specified this calculation 
is by-passed. 


The acentric factor w is computed from Edmister's equation, 


log (P./14.7) 


3 
vo => ——— - 1.0 (4-26) 
i Tiles =i) 0 
where T is mean average boiling point in °R. 


b 


The latent heat of vaporization at the normal boiling point 


is calculated by Kistiakowsky's equation, 


= = 7.58 + 4.571 log Ty (4-27) 


and adjusted to 25°C by means of Watson's equation, 


AH. T_ - 537 0.38 
re25Rby lg agEc ) (4-28) 


AH Ps - qT 


The solubility parameter is calculated by, 


(AH - RT)? 
R ee _penress (0-3 eee tees (4-29) 
| v_ 


The constants required for the ideal gas state enthalpy 


equation, 


He = (AH) .+(BH) T+ (CH) T+ (DH) ,T° (4-30) 
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are calculated from the following equations, 


os 2 3 2 3 
AH MW(Co+C A+CA +04 #0), KCK +k +C  gAK) (4-30-a) 
BH = MW(C_+CoA+C Natt A+ KC K-4¢ Kk?) (4-30-b) 
dm oF £9 10 1] 12 13 
as . 
CH = MW(C), + C1 5A) (4-30-c) 
DH = MW(C, ¢ + CJA) (4-30-d) 
where 
A = API 
K = the UlO2P. characterization factor 


MW = Molecular Weight. 


4.4 Limitations of the Chao-Seader Correlation 
The general Chao-Seader correlation has been widely used and 
extensively evaluated. The method is characterized by a generalization 
based on sound thermodynamic relations. Furthermore the non-ideality 
is related to the properties of pure components only. Then the 
calculation procedures are greatly simplified. On the other hand the 
existence of some limitations due to the simplified assumptions is 
indispensable. Lenoir (46) has summarized the conditions under which 
the correlation can give satisfactory predictions as follows, 
(1) The pressure should not exceed 1000 psia and should not 
be greater than 0.8 of the system critical pressure. 
(2) The reduced temperature of any hydrocarbon other than 
methane should be between 0.5 and 1.3. 
(3) The methane content in the liquid phase should not exceed 
0.3 mole fraction, and for any system containing hydrogen 
or methane, the reduced temperature should not exceed 0.93, 


and the temperature should be between 100 and 500°F. 
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(4) When predicting paraffins or olefinic equilibrium ratios, the 


aromatic content should be less than 0.5 mole fraction. 


4.5 Partial Derivatives of K-ratio and Enthalpy 
The general algorithm for multistage multicomponent separation 
calculations proposed by the author requires the following partial 
derivatives, 
OK. j OK, j oH oh 
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Either numerical or analytical techniques may be used to 
evaluate these partial derivatives. The computer effort is generally 
large for numerical evaluation. As far as complicated correlations 
such as the Chac-Seader method are concerned, it is especially 
advantageous to determine partial derivatives analytically. 


K-ratio is expressed in the form, 


Then the partial derivatives of K with respect to T and 
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where analytical expressions for the derivatives of the vapor phase 
fugacity coefficient, Vis the liquid phase fugacity coefficient, 
/P, and of the activity coefficient, Y;» can be derived from 


the relations making up the Chao-Seader correlation. 
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Now define f(T) and g(T) as follows, 
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temperature and pressure. 


The partial derivatives of the activity coefficient are 


evaluated as follows, 
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However, SES) eT ee and repay ees be taken to be 


negligible, since constants A. , A, B. and B are functions of vapor 


phase compositions, y;'s and those constants are approximated to 
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is expressed in terms of partial derivatives of the constants with 


respect to x, 'S. Consequently, 
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CHAPTER 5 


THE NEW METHOD 


A new general calculational method for equilibrium stage 
processes is described in this chapter. The method employs the Taylor 
approximation procedure for the total linearization of all the model 
equations. This linearization technique, which is actually a modified 
Newton-Raphson method, has been widely applied for solving nonlinear 
equations numerically. For example, Lee (47) has solved successfully 
various kinds of nonlinear differential equations. The present 
author (48) applied the technique to solving complicated partial 
differential equations which represent the reactor problem where 
chemical reaction takes place between gas and liquid reactants on 
solid catalysts. In this case simultaneous heat and mass transfer with 
fluid mixing was taken into account. Roche and Staffin (49), (50) 
applied this linearization technique to calculations for a ternary 
liquid - liquid extraction problem. The present author (51) extended 
the same concept to a multicomponent extraction problem. Recently 
Gentry (7) has proposed the application of this linearization 
technique to distillation processes. 

In the method proposed by the author all the model 
equations, namely, over-all material balances, component material 
balances, energy balances and the summations of compositions, are 
approximated in a linearized form. The linearized equations are 
then solved simultaneously, without decoupling, for corrections to 
temperatures, flow rates and compositions. In Tomich's method only 


the energy balance and summation equations are linearized decoupling 
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the material balance equations. The methods are significantly 
different in that the linearized equations of the method proposed 

by the author are expressed analytically. However, those of Tomich's 
method cannot be expressed analytically since the functional relation- 


ships of the equations with related variables are implicit. 


5.1 Linearization of Equations 
A multivariable nonlinear function, f, is often approximated 


as follows, 


f(x) +Ax,, X,tAX, a x tAx) 
n 
of 
~ F(K) X59 eX) ot +), ae : (5-1) 


Equations (2-3), (2-4), (2-5) and (2-6) are approximated 


in the same manner as equation (5-1), then linearized setting, 


f(x) +Ax,, X5tAXs » 


ae x thx.) =O . (5-2) 
However, slight modification of standard linearization 

procedures will be utilized in this thesis to obtain convenient 

matrix structures and manipulation. Namely, only the partial 

derivatives that have a dominant influence on the solution are included 

in the linearized equations. The details will be discussed later. 


The linearized equations are: 


Over-all material balances 
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ie ite j j j : ( 
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AL = } My + AV we Gal 2>N (5-4) 
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AN ae aL VAY (= a (5-5) 


is specified 


and AV., = AL, = AV, = Mie 0 for a distillation problem when V 


where AV, = 0 for a reboiled absorber problem when V 


l 


and Vo are specified. 


ee ee eee (5-6) 


where ALY = 0 except for the usual absorber problem. 


Component material balances 
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AQ, = 0 for an absorber problem 
AV, = AQ, = 0 for a reboiled absorber problem 
AV) = AL, = AV. = Q for a distillation problem. 
At j =N 
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N-] N N 
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+ ALY hy )-AV\Hy"ALyhy AO, = -Ey (5-12) 


where 


45 


> 
/=) 
i} 


0 for an absorber problem 


> 
i 
il 


N 0 except for an absorber problem. 


Summation equations 


yiGAx.OF sneSyt} jschllsaNs (5-13) 


The equilibrium ratios, and the liquid and vapor enthalpies 
are functions of temperature, pressure and phase compositions and may 
be calculated as such at each iteration. However, for the purpose of 


linearization only, the following relations were assumed. 
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The number of partial derivatives to be evaluated follows. 
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Then the total number of partial derivatives is 2N x (1+NCP). 

The linearized equations then only include the partial 
derivatives that have the strongest influence on the solution and this 
in turn greatly reduces the computer storage required. If the standard 
linearization procedure is employed directly without modification, 
simple matrix structures as in this thesis cannot be obtained and the 


matrix manipulation will encounter prohibitive difficulty. However, 
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this practical simplification of the linearized equation need not 
affect the rigor of the final result in that the equilibrium ratios 
and enthalpies are calculated at each iteration as functions of 
temperature, pressure and normalized phase compositions. 

Either numerical or analytical techniques may be used to 
evaluate the partial derivatives which appear in the linearized 
equations. The amounts of computer storage and calculational time 
required are generally large for numerical evaluation. The computer 
program developed in this thesis employs the Chao-Seader correlation. 
The partial derivatives are evaluated analytically in the program. 
The analytical expressions for the partial derivatives are given in 
Chapter 4. 

Equation (5-4), the linearized over-all material balance, 
may be substituted into the component material balance equations 
and the energy balance equations to eliminate all as terms. The 
linearized component and energy balances may then be expressed in 


matrix form as follows. 


Component material balances for an absorber problem 
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Component material balances for a reboiled absorber where V is 
specified. 


Equation (5-17) is modified as, 
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AV, = 0 


Component material balances for distillation where V, and V, are 
specified. 


Equation (5-17) is modified as 
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Energy balances for an absorber 
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Energy balances for a reboiled absorber 
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where p,q, r, S, t, and u are the same as those in equation 
(5-18-a). 
Energy balances for distillation 
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where p,q, r, S, t and u are the same as those in equation 


(5-18-a). 
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Now the component material balances combined with the 
equilibrium relationships and the energy balances are expressed 


simultaneously in matrix forms as equations (5-17) and (5-18). 


5.2 cOlution Procedure 

The solution procedure is iterative. Figure 2 shows the 
flow diagram for computation. As seen in the figure all unknown 
variables of temperatures, flow rates and phase compositions are 
solved simultaneously in a single convergence loop. This is a 
significant characteristic of the new method. 

At first appropriate specifications are made and initial 
values are selected for the flow rates, aie and V,'s, the 
temperatures, ie and the phase compositions, x. .'s_ and 


VyJ 


Y; °° Substitution of ideal K-ratios into equation (2-6) and 


solving the C-matrix with C. : = 0 by the Thomas algorithm such 


> 


as described by Wang and Henke (4) will readily provide initial 


estimate of xX; ek Then the phase compositions are normalized 
3K oK 3H oh : 
STatel> ye a ais OT? Bx? BT? and sy are calculated as described 


in Chapter 4. All the elements in the matrices of equations (5-17) 
and (5-18) can then be calculated. From equations (5-17), (5-18) 
and the summation equation (5-13) with oF = 0, the correction 
terms for the unknown variables are obtained. If the corrected 
Variables satisfy a cegtain convergente criterion, then the 
computation will stop. Otherwise the procedure will be performed 
again iteratively. 

The detail computation procedure for obtaining the solution 


is as follows, 
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PL GURE 2 


FLOW _ DIAGRAM FOR COMPUTATIONAL 
PROCEDURE 
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and 


where 
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Greene ila 
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Convert the first matrices in equations (5-17) and (5-18) 


to unit matrices and obtain matrix equations of the form, 


[ei + Blt CAVE one (5-19) 


| AT+ PAV=Q . (5-20) 


The following algorithm was derived for inverting 


the tridiagonal matrices (51). 
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However, one could also use the Gauss-Jordan elimination 
procedure. The algorithm will be simplified by the use of 
the new procedure. A numerical example which illustrates 
the application of this procedure is shown in Appendix 2. 


(2) Sum equation (5-19) over i = 1 > NCP 


1ZAx. + BAT +CAV=D . (6-29) 


(3) Eliminate the Ax, term from equation (5-19) by substituting 
equation (5-13) with S$. =0 to obtain, 


==) eee (5-23) 


(4) Substitute equation (5-20) into equation (5-23) to obtain 
ee and then Mie are calculated from equation (5-20). 
The values are substituted into equation (5-17) which is 


then solved for the Ax, ve using the Thomas algorithm. 


> 


(5) Compute new Vi's, TyIs and x, ;° ase 
een sen (5-24) 
J J J 
iageaclr—saa aie ais (S55) 
a | J J 
Ah = os . PATA (5-26) 
154} I,J I,J 
where t, (0< t <1) is a weighting factor. 
A weighting factor, t, is chosen so as to satisfy 
the following relation. 
k+] k 
[CRIT] < [CRIT] (5227) 


where [CRIT] is defined as, 
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(6) Y; j 3 are obtained from equation (2-1) and then normalize 
> 


the phase compositions, x; i and Y; gee Theoretical ly 


> 


the normalization of x, .'s is not required. However they 


I, 


are normalized to avoid truncation error. Again calculate 
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These values and the new stage variables are used to 
obtain the elements in the matrices. 
(7) Continue calculations until a convergence criterion such 


as the following is satisfied. 


ee (5-29) 


One possible definition of the convergence criterion may be, 


N NCP 
ce) CoA Be aecte & of ; (5-30) 


ey ; 


This represents the summation of imbalances in over-all material 
balances (2-5), component material balances (2-6), energy blances 
(2-3) and summation equations (2-4) over an entire column. However, 
in the computer program flow rates are calculated so as to satisfy 
equation (2-5) with M; = 0. In addition, the liquid phase 


.'s always satisfy equation (2-4) with Sj S70 


compositions, x; ; 
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Then equation (5-30) is expressed as, 
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In place of equation (5-31) the convergence criterion may 


be defined as, 
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NCP 


N 
ferries mel tec Shey, for pAanh pete (5-32) 
ia alte de : 


Since If the following conditions are satisfied, 


then the component matrerial balance equations should be satisfied, 
namely, ey = 0 

The second term in equation (5-32) is normalized by dividing 
by the total heat input to any stage j so that it is of the order 
of magnitude 1. Consequently equation (5-28) can be defined as a 
convergence criterion. 

The weighting factor, t, in step (5) is chosen as follows. 


chy 


First a weighting factor, at each iteration is determined so 


as to satisfy the following conditions. 
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0.5 vi a “ taV <2 0 ve ; j lene (5-34) 
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Then .t is the smallest value of t which satisfies the above 
conditions. If the condition (5-27) is not satisfied, the next factor 
is Set as (2) = 0.5 £6), lt has been the author's experience that 
the condition (5-27) is usually satisfied with the first weighting 
factor. p61), The conditons (5-33) and (5-34) are selected so as 

to give a steady solution. For example, in case the temperature 


correction at any stage is over 1000°F, obviously one cannot expect 


a converged solution. 
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some flexibility in speerfications extsts in’ that the 
algorithm for an absorber may be used for distillation problems where 
all heat loads including the condenser and the reboiler are 
specified, and for a reboiled absorber where the reboiler load is 


speci fied. 


CHAPTER 6 


EVALUATION OF THE NEW METHOD 


The calculational procedure described in Chapter 5 has been 
used to prepare a general program for the solution of multicomponent, 
multistage separation problems. The new algorithm has been tested on 
a number of problems including absorbers, reboiled absorbers and 
distillation columns typical of those used in natural gas processing 
industries, and complex hydrogen rich demethanizers used at the 
Canadian Industry Limited polyethylene piant in Edmonton and at the 
Polymer Corporation Limited synthetic rubber plant in Sarnia, Ontario. 
Comparisons have been made with the results obtained using algorithms 
prepared by the author based on Tomich's method (6) and the 
tridiagonal matrix-Bubble Point procedure described by Wang and 
Henke (4), and with the results obtained by Burningham and Otto (22) 
for programs based on (i) the Sum-Rates-Thiele-Geddes procedure, 
(ii) the Sum-Rates-Tridiagonal matrix procedure and (iii) the 
Bubble Point-Thiele-Geddes method with 6 convergence. 

Composition dependent equilibrium and enthalpy data were 
calculated directly in all programs using the Chao-Seader correlation. 
Computations were made on the |.B.M. 360-67 system at the University 
of Alberta. 

Descriptions of the problems utilized in evaluating the 
new algorithm follow. Complete computed results obtained with the 


new algorithm are given in Appendix 3. 
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6.1 Test Problems 


6.1.1 Example Problem No. 1 - Absorber 
This absorber problem is described by Holland et al (2). 


The problem specifications and initial assumptions are as 


follows. 


Compositions of Entering Streams 


Component Feed Lean Oi] 
(moles/time) (moles/time) 
methane 70 0 
ethane AS 0 
propane 10 0 
n-butane ky 0 
n-pentane | 0 
n-octane 0 20 
Pressure = 300 psia 
Number of theoretical stages = 8 
Lean oi] temperature = 90° F 


Feed temperature = dew point 


Initial Assumptions 


Stage No. Temperature Vapor Rate 
(oF) (moles/time) 

| | 100.0 85.0 

2 100.0 9020 

3 100.0 90-0 

h 100.0 90.0 

5 100.0 90.20 

6 100.0 90@0: 

7 100.0 9050 

8 100.0 95.0 

6.1.2 Example Problem No. 2 - Absorber 


This problem is described by Burningham and Otto (22), and 
is typical of absorbers in natural gas processing plants. 
The problem specifications and initial assumptions are 


as follows. 
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Compositions of Entering Streams 


Component Feed Lean Oi] 
(moles/time) (moles/time) 

nitrogen 20627 0.0 
CO, 1s235 O59 
methane 8721.5 0.0 
ethane IWAey fats! 0.0 
propane 686.8 0.0 
i-butane 12.8 0.0 
n-butane 136.7 0.0 
i-pentane 19.7 0.0 
n-pentane 19.5 0.0 
n-hexane [33 0.0 
hypth-1 8.4 027 
hypth-2 | 2.5 229 
hypth-3 0.5 129.9 
hyp th-4 0.1 204.3 
hypth-5 0.1 guy 

Pressure = 1310.0 psia 

Number of theoretical stages = 6 

Lean oi] temperature = -5.0° F 

Feed temperature = 45.0° F 


Inter cooler duty on plate 6 = ~-5,800,000.0 BTU 


alt 


Hypothetic 
Normal B 
Point ¥{ 
hypth-] 200 
hypth-2 265 
hypth-3 325 
hypth-4 380 
hypth-5 425 
Initia 
Stage Temp 

| 

ri 

3 

h 

5 
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Example Problem No. 3 - 


Reboiled Absorber 


al Component Data 

oiling AP | Molecular 

=e) Gravity Weight 
6352 89.6 
B73 103.6 
52.8 12720 
47.3 145.0 
Lh 6 166.0 

] Assumptions 

erature Vapor Rate 

(2F) (moles/time) 

2070 9398.6 

290 9500.0 

25.0 72.0 

25.0 10000.0 

25.0 10500.0 

20.0 11000.0 


This problem is described by Holland and et al a) 


The specifications and initial assumptions are: 
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Compositions of Entering Streams 


Component Feed Lean Oi] 
(moles/time) (moles/time) 
methane 80.0 0.0 
ethane 6.7. 0.0 
propane 6.7 0.0 
n-butane 6.7 ; 0.0 
n-octane 0.0 30.0 
Pressure = 300 psia 
Top product (vapor) = 95.0 moles 
Liquid side stream from sie plate = 15.0 moles 
Lean oil temperature = 90°F 


Feed temperature = dew point 
Number of theoretical stages = 10 


Feed stage number = 6. 


Initial Assumptions 


Stage No. Temperature (°F) Vapor Rate 
(moles/time) 


100.0 95.0 
2 150.0 105.0 
3 150.0 105.0 
5 150.0 110.0 
6 220.0 75.0 
7 280.0 50.0 
8 350.0 50.0 
gs) 400.0 50.0 
10 480.0 65.0 


rnin w= 


_ 7 e- 


6.1.4 Example Problem No. 4 - Reboiled Absorber 

This reboiled absorber problem is described by Burningham 
and Otto (22), and is typical of reboiled absorbers utilized in 
natural gas processing plants. 


The specifications and initial assumptions are: 


Compositions of Entering Streams 


Component Feed Lean Oi] 
(moles/t ime) (moles/time) 

nitrogen 0:7 0.0 
CO, 12.4 0.0 
methane 167.4 0.0 
ethane 474 LS 0.0 
propane 4hLO.2 0.0 
i-butane 64.3 0.0 
n-butane 128.8 0.0 
i-pentane 19.5 0.0 
n-pentane 19.4 0.0 
n-hexane leu 0.0 
hypth-1 9.0 US) 
hypth-2 ae: 7A 
hypth-3 20005 SES} 
hypth-4 20373 esses) 
hypth-5 24.7 Tent 

Pressure = 270.0 psia 

Number of theoretical stages = 16 

Lean oi] temperature = -5.0°F 

Feed temperature = 50.0°F 


Side reboiler duty on plate 12 = 4,000,000 BTU 
Feed stage number = 6 


Top product (vapor) = 677.0 moles 


Hypothetical Component Data 
OO ete neha ASB eS Sh 


Normal Boiling API Molecular 

Pointe (oF) Gravity Weight 
hypth-1 200 63.2 89.6 
hypth-2 265 5723 103.6 
hypth-3 325 52.8 12720 
hypth-4 380 47.3 145.0 
hypth-5 425 4h 6 166.0 


impel Assumptions 


Stage No. Temperature Vapor Rate 

me ce (moles/time) 
l ene 677.9 
2 90.0 750.0 
3 100.0 800.0 
4 110.0 850.0 
5 120750 900.0 
6 135..0 3950.0 
H 15020 900.0 
8 160.0 900.0 
9 175.0 900.0 
10 190.0 900.0 
1] 200.0 9507.0 
12 210.0 1000.0 
13 22020 1000.0 
14 230.0 1000.0 
15 240.0 1000.0 
16 250.0 1000.0 
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6.1.5 Example Problem No. 5 -  Reboiled Absorber 

This is the hydrogen rich demethanizer operating at the 
Canadian Industry Limited polyethylene plant in Edmonton. 

The problem specifications and initial assumptions are 


as follows. 


Compositions of Entering Feed 


- Component Feed Lean Oi] 
(moles/time) (moles/time) 

methane 12777920 0.0 
ethane 278.870 0.0 
ethylene 413.190 0.0 
propane VIZ/9 0.0 
hydrogen 0.0 0.0 
benzene 0.0 792.610 
toluene 0.0 70.980 
xylene 0.0 23.660 
propylene 0.0 592150 
ethylcyclohexane 0.0 23.660 

Pressure = 464.7 psia 

Number of theoretical stages = 20 

Lean oil temperature = 5°F 

Feed temperature = 50°F 


Inter cooler duties 


on plate 4 = 840,000.0 BTU 
on plate 7 = 1,095,000.0 BTU 
Feed stage number = 17 


Top product (vapor) = 642.2 moles 


Initial Assumptions 


Stage No. Temperature Vapor Rate 
ees SS (°F) (moles/time) 
l 0.0 642.0 
a 5.0 650.0 
3 10.0 650.0 
mn 5.0 700.0 
5) 10.0 900.0 
6 15.0 920.0 
7 10.0 930.0 
8 1220 1100.0 
9 20.0 15020 
10 30.0 1200.0 
i] 40.0 1200.0 
12 5050 1250.0 
13 60.0 1250.0 
14 70.0 1300.0 
15 715.20 1300.0 
16 80.0 1350.0 
17 90.0 1400.0 
18 120.0 Zo0. 
ibe 130.0 350.0 
20 15050 400.0 
6.1.6 Example Problem No. 6 - Reboiled Absorber 


This problem is hydrogen rich three feed demethanizer 
Operating at the Polymer Corporation Limited synthetic rubber plant 


in Sarnia, Ontario. It is taken from the information given by Prowse 


a re 7. a = S eae = a ee | 


= 
— 


ee ea 
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and Johnson (27), and Petryshuk and Johnson (52). 
The problem specifications and initial assumptions are as 


follows. 


Compositions of Entering Streams. 


Component Feed - | Feed - 2 Feed - 3 
(moles/time) (moles/t ime) (moles/t ime) 
hydrogen 0.0 0.0 562.4 
nitrogen 0.0 0.0 207.2 
methane 0.0 28.9 S| ar 
ethylene 0.0 Bye yh. 
ethane 0.0 144.5 Kya 4 
propylene 0.0 491.3 384.8 
propane 23.66 | 317.9 11324 
i-butene 196.56 520.2 88.8 
2-butene 4h 08 173.4 14.8 
n-butane Vi bee ta: 982.6 Se Ie 
n-pentane 331.24 158.95 35.52 
heptane 47.32 14.45 8.88 


Pressure = 475.0 psia 


Number of theoretical stages = 30 
Lean oi] temperature = -4°F 

Feed - 1 temperature = -4°F 

Feed - 2 temperature = -4°F 

Feed - 1 stage number = 15 

Feed - 2 stage number = 22 


Top product (vapor) = 1520.0 


The initial temperatures are assumed to be linear between 


O°F and the top stage and 203°F at the reboiler. 


Initial Assumptions 


Stage No. Temperature Vapor Rate 
l 0 1520 
z 7 1600 
3 14 1650 
4 21 1700 
5 28 1750 
6 35 1800 
i 42 1850 
8 kg 1900 
9 56 1950 

10 63 2000 
11 70 2000 
t2 yer 2000 
13 84 2000 
14 91 2000 
15 98 2000 
16 105 2000 
7, Hike 2000 
18 yale! 2100 
1g 126 2200 
20 133 2300 
21 140 2400 
o¥8) 147 2500 
23 154 1000 
24 161 1200 
25 168 1400 
26 liz 1600 
27 182 1800 
28 189 2000 
29 196 2200 


203 2400 
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6.1.7 Example Problem No. 7 -. Distillation Column 


This distillation column is described by Wang and Henke (4). 


The problem specifications and initial assumptions are: 


Compositions of Entering Streams 


Component Feed - | Feed - 2 
(moles/t ime) (moles/time) 
methane 2.0 0.0 
ethane 1O70 0.0 
propylene 6.0 1.0 
propane 12.0 720 
i-butane 30 4.0 
n-butane 3.0 1720 
n-pentane O75 V5e2 
n-hexane 0.0 920 
heptane 0.0 4.5 
octane 0.0 4.3 
decane 0.0 3m 


Pressure = 264.7 psia 
Number of theoretical stages = 2] 
Feed stage numbers = 7 and 13 


Feed temperatures 


stage 7 = dew point 
stage 13 = bubble point 
Vapor distillate rate = 23.0 moles 


Reflux ratio 22 3.122 


Liquid side stream from yen plate 15.0 


2520 


th 
Vapor side stream from 16 plate 


ae 


-s 


a 


= 


ritual 


Initial Assumptions 
ee ELON 


Plate No. Temperature 

(SF) 

I 80.0 
2 95.5 
3 70 
4 135.5 
5 154.0 
6 17225 
7 191.0 
8 209.5 
2 228.0 
10 246.5 
1 265.0 
12 265.0 
13 302.5 
14 320.5 
15 339.0 
16 Bey ie: 
1, 376.0 
18 394.5 
19 413.0 
20 431.5 
21 450.0 


Vapor Rate 
(moles/time) 


23. 
94. 
4. 
. 806 


94 


Qh. 


Qh 


94. 
60. 
60. 
60. 
60. 


60 


60. 
60. 
60. 
60. 
85. 


85 


85. 
85. 
85. 


Gale Example Problem No. 8. Distillation Column 


This is a two product condenser stabilizer which is presently 


0 
806 
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306 
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Operating in a natural gas processing plant in the province of Alberta. 
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The problem specifications and initial assumptions follow. 


Composition of Entering Streams 


Component Feed 
(moles/time) 


CO. 0.419 
nitrogen OEZ2ul 
methane 74.737 
ethane 78.890 
propane 130.060 
i-butane Ah 244 
n-butane 90.173 
i-pentane 41.344 
n- pentane 43.606 
2-methyl pentane 15.143 
hexane 14.955 
heptane 25.017 
octane 15.875 
nonane 8.784 
decane Zeoge 
undecane 0.443 
dodecane 0.174 
tridecane 0.088 


Pressure = 350.0 psia 

Number of theoretical stages = 25 
Feed stage number = 9 

Feed temperature = 250.0°F 

Vapor distillate rate = 21.02 moles 
Liquid distillate rate = 21.02 moles 


Reflux ratio = 1.9 
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Initial Assumptions 


Stage No. Temperature Vapor Rate 
| 120.00 210.20 
2 131267 1219.20 
3 143.75 1219.20 
4 155.62 1219.20 
5 167.50 1219520 
6 ey) 1219.20 
7 191.25 1219.20 
8 20312 1219.20 
) 215.00 1219.20 

10 226.6] 800.19 
1] 235.95 800.19 
12 250.62 800.19 
13 263.40 800.19 
14 274.37 800.19 
15 286.25 800.19 
16 298.12 800.19 
LF/ 310.0 800.19 
18 EPAV Tw / 800.19 
19 333.d0 800.19 
20 345.62 800.19 
21 357.50 800.19 
02 369.37 800.19 
rae 381).25 800.19 
2h 393.12 800.19 


25 405.00 800.19 
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6.1.9 Example Prob] 
This is a 


province of Alberta. 


em No. 9 - Distillation Column 


n-butane-i-butane splitter operating in the 


A linear temperature profile and constant molar 


flow rates are assumed as initial assumptions. 


The problem specifications are: 


Feed Composition 


Component Feed 
(moles/time) 


propane 10.0 
i-butane 652.0 
n-butane 1140.0 
i-pentane 1735 
Pressure = 84.0 psia 
Number of theoretical stages = 70 
Feed stage number = 35 
Feed temperature = boiling point 
Vapor distillate = 670.0 


Reflux ratio = 11.3 


G22 Comparison of the Results 


The computed results obtained with the new algorithm are 


shown in Appendix 3. 

required are summari 
In Table | 
New: 
Tomich: 


BE PR pe: 


The computing time and number of iterations 
zed in Table | and Table 2. 
and Table 2, the methods compared are, 
New method proposed by the author. 
Tomich's method. 
Bubble Point method combined with the Tridiagonal 


Matrix procedure. 
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SRG: Sum Rates method combined with the Thiele- 
Geddes procedure. 

Shei Ri: Sum Rates method combined with the Tridiagonal 
Matrix procedure. 

BPS TG. Bubble Point method combined with the Thiele- 


Geddes and the Theta convergence procedure. 


In Table 3, Table 4 and Table 5, typical comparisons of the 
rate of convergence between the new method and the Tomich method are 


shown. 
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Table | 


Comparison of Computing Time 


(Time is seconds on an IBM 360-67) 


Absorber 


(79) 


A 


14.43 


IM ES 


Reboiled Absorber 


86.72 


pe De ee 19009 | S431 08 


| 
g | 21.35 | 123.84 


Distillation 


‘ng ie 


AA: 


Table 1 Comments 


divergent 

on an IBM 7040 

Burningham and Otto (22) 

Otto, private communication 

not converged in specified time (180 sec.) or number 


of iterations (30 times) 


Wang and Henke (4); they reported that BP-TRI is 
faster than BP-TG for problem 7 when polynomial 


data are employed. 
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Table 2 


Comparison of Number of Iterations Required 


TOMICH | BP-TRI SR-TG | SR-TRI BP-TG 
problem 


Absorber : 
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X: divergent 
not converged in specified time (180 sec.) or number 
of iterations (30 times) 


see the comment in Table | 
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Table 3 


Comparison of Rate of Convergence with 


Tomich's Method 


(Problem 1) 

The New Method The Tomich Method 
Iteration [CRIT] Iteration [CRITI 
Sono eee ___ No. 

0 9.324% 10°! 0 0.453 x 107! 
0.361 x 10> ovetgex 107+ 
2 0.145 x 10°° 5 0.633 x 107° 
=5 = 

3 On7o2 x 10 10 0.583 x 10 
a! 

ie 0.343 x 10 


20 0.908 x 10 
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Table 4 


Comparison of Rate of Convergence with 


Tomich's Method 


(Problem 3) 


The New Method The Tomich Method 
Iteration TCRIT] Iteration [CRIT] 
eee Nove iter Noten 

0 0.545 0 0.424 

1 0.197 0.270 

2 Deere 100. 5 0.119 

3 0.122 x 1077 10 0.434 x 10 

4 voce on 15 0.289% 105 

5 eye ee 20 15505" 10 
25 0.139 x 10 


24 0.145 x 10 
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Table 5 


Comparison of Rate of Convergence with 
a 


Tomich's Method 


(Problem 7) 
The New Method The Tomich Method 
Iteration [CRIT] Iteration [CRIT] 
No. No. 
0 0.873 0 0.983 
0.199 0.675 
2 0.198 x 1077 5 One 2 ano 
3 0.461 x 10°? 10 0.106 x 10> 
I i 


h 0.140 x 10- 14 0.168 x 10. 


6.3 Discussion 

A comparison of the performance of the new method with 
Tomich's method and with Bubble Point and Sum Rates procedures is 
given in Table 1 and Table 2. The Sum Rates and Bubble Point 
procedures are based on traditional procedures and, in general, are 
not reliable for all types of problems. The Sum Rates procedures 
are more effective for absorber problems whereas the Bubble Point 
methods are more reliable for distillation problems. These 
observations are in agreement with the conclusions of Friday and 
Smith (14). 

The Thiele-Geddes bubble point method with theta convergence 
has successfully solved some complex reboiled absorber problems. 
However there is no assurance that this method will converge for a 
particular problem (52), (22). 

Both Tomich's method and the new method show promise of 
being reliable for all types of separation problems. However, 
Tomich's method requires longer computational time and has 
difficulty with complex reboiled absorber problems. 

Demethanizer absorbers which handle a significant amount 
of hydrogen, such as described in problems 5 and 6, seem to be one 
of the more difficult types of problems to solve. Petryschuk and 
Johnson (52) were not able to solve a problem similar to example 6 
with their program based on the Thiele-Geddes bubble point procedure. 
Prowse and Johnson (27) solved a problem similar to example 6 using 
a method based on Ball's relaxation procedure (25). Relaxation 
methods are characteristically stable but converge quite slowly. 


Prowse and Johnson (27) report using more thaneG0 itera ttons ) forta 
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problem similar to example 6. 

Problem 5 was formulated from information provided by C.1I.L. 
(Canadian Industry Limited) for a column Operatvingeaternecet sini. 
Edmonton polyethylene plant. The new method readily solves this 
complex column. However, attempts to simulate the column operation 
have, as yet, been unsuccessful because of uncertainties as to the 
identity of all the components present, the plant analyses and as to 
the proper equilibrium and enthalpy data. 

The results obtained for problems 7, 8 and 9 indicate that 
the new method is the superior method of the procedures tested for 
distillation problems. The method was particularly effective for the 
iso-butane -n-butane splitter which contains a large number of 
equilibrium stages and is described in problem 9. 

Comparisons of the convergence characteristics of the new 
method and of the Tomich method are shown in Table 3, Table 4 and 
Table 5. In general, the new method requires less computational 
time and fewer iterations than the Tomich method. However, the 
average computational time per iteration is less for the Tomich 
method. This illustrates the effectiveness of the Broyden 
procedure. 

The initial assumptions of temperature and flow rate 
profiles can have a marked effect on the performance of any procedure, 
thus identical initial assumptions must be used for valid comparisons. 
The initial assumptions used for the test problems were moderately 
far from the final solutions and were selected as ones that could be 


easily calculated or assumed without prior knowledge of the problem 
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solution. The new method shows promise of being relatively insensitive 
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to the initial assumptions used, that is, it appears to give convergent 
solutions for a wider range of initial assumptions for any given 
problem than other methods. However, proper analysis of the relative 
sensitivity of the methods to initial assumptions is complex and must 
take into account the method of finding the weighting factor t. 
Additional work is required to clarify the effect that the choice of 
initial assumptions has on the convergence characteristics. 

The method of selecting the weighting factor, t, has an 
important effect on the characteristics of the new method and Tomich's 
method. Further work is required to provide the optimum procedure 
for selecting the weighting facotr. 

The new method obviously requires a larger amount of computer 
storage than the traditional decoupled methods such as those 
proposed by Lewis and Matheson (1), Thiele and Geddes (3), and Wang 
and Henke (4). However, singinficantly less computer storage is 
required than for the methods proposed by Gentry (7) and Naphtali et al 
(9). The methods of Gentry and Naphtali et al require a prohibitive 
amount of computer storage when the use of composition dependent 
equilibrium and enthalpy data is incorporated in the methods. The 
new method requires less storage than the Tomich method. A 
comparison of the sizes of the inverses of the Jacobian matrices can 
be used as a measure of the computer storage requirement. For example, - 
for a 100 stage column the total number of elements in the inverse of 
the Jacobian matrix for the new method is 10" whereas 4 x 10" elements 
are required in the Tomich method. 

The Chao-Seader correlation has been used exclusively for 


the direct calculation of equilibrium and enthalpy data for the 
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problems reported in this thesis. Other correlations such as the ones 
proposed by Chueh and Prausnitz (38), Lee and Edmister (40), Wilson 
(33) as well as polynomial data correlations can be easily used with 
the new method. 

It should be possible to develop analytical expressions for 
the partial derivatives required from the correlations mentioned above. 
Numerical evaluation of the partial derivatives will be required when 
analytical expressions are not available. 

When linearizing the model equations, only the partial 
derivatives that were thought to have a dominant influence on the 
solution were employed. This concept is fundamental to the new 
method in that it allows convenient matrix structures to be formulated 
which in turn facilitate the simultaneous solution of all the model 
equations. The test of the method on a wide range of problems 


suggests that this strategy is quite successful. 
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CHAPTER 7 


CONCLUSIONS 


The new calculational procedure for multistage, multicomponent 
separation calculations that has been described and evaluated in this 
thesis is a reliable and efficient method for absorber, reboiled 
absorber and distillation calcualtions typical of those encountered 
in the hydrocarbon processing industries and the chemical industries 
in general. 

The new method is more reliable and requires less com- 
putational time than the Tomich method and Bubble Point and Sum Rates 
procedures, and would appear to be a significant improvement on 
current procedures. 

The method seems to be relatively insensitive to initial 
assumptions, but additional work is required to clarify the influence 
of initial assumption and of the choice of the weighting factor on 
the convergence characteristics. 

Structuring the equilibrium stage problem as described in 
this thesis, making effective use of matrix manipulation procedures 
to solve all the model equations simultantously, avoiding iterative 
bubble point calculations or other decoupling methods and the fact 
that only 2N(I+NCP) partial derivatives need to be evaluated per 
iteration, all contribute to the significant reduction in computer 


storage requirement and computational time compared to other methods. 
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AH, BH, CH, DH 


arb, c 
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xi 


AH 


NOMENCLATURE 


API gravity, or Chao-Seader coefficient 
defined by equation (4-3) 

Jacobian matrix | 

ideal gas constants 

elements of matrix 

matrix 

function defined by equation (2-2) 
matrix 

matrix 

densi ty 

elements of matrix 

function defined by equation (2-3) 
normalized energy imbalance 

general function 

fugaci ty 

defined by equations (4-35) and (4-36) 
elements of matrix 

column vector defined by equation (3-13) 
molar feed rate 

element of matrix 

enthalpy of liquid, or defined by equation (4-2) 
enthalpy of vapor 

matrix defined by equation (3-19) 

heat of vaporization 


unit matrix 
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equilibrium K-ratio or U.0.P. characteristic factor 


molar liquid flow rate 

function defined by equation (2-5) 
molecular weight 

element of column vector 

total number of stages 

number of components 

column vector defined by equation (3-15) 
elements of matrix 

pressure 

matrix 

heat duty 

column vector 

gas constant 

elements of matrix 

function defined by equation (2-4) 
molar flow rate of liquid side stream 
molar flow rate of vapor side stream 
weighting factor 

temperature 

column vector 

element of column vector 


molecular volume of gas 


molar vapor flow rate, or liquid molar volume 


column vector 
liquid composi tion 


colum vector defined by equation (3-13) 


(21-€) noltaaps yt bShNASb) wae". omstoo 
sweesrg 

atten 

yaub stood 

--9etI0N ambos 

; Jewsenos exp 

aidan to adnensis 

(0-2) nolseupa yd banlieh notysnut 


messse sbie biugpl! lo eter wot? isle 
etrita: eb Fe 20qey Necaney eet? vetos 
| sais? ar)trigisw 


~~ 


©. 


magis,5 £554 )3 « 


column vector whose element is £ Ax 


column vector 
vapor composition 
column vector defined by equation (3-18) 


feed composition, or compressibility factor 
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Greek Letters 


o element of matrix in equation (5-17) 
8 element of matrix in equation (5-18) 
Y activity coefficient 
6 solubility parameter 
A "new vatue minus old value!! 
€ defined by equation (5-21) 
v liquid fugacity coefficient 
o Euclidean norm 
Vy) vapor phase fugacity coefficient 
mn) acentric factor 
Subscripts 
F feed 


i component number 


j stage number 
L liquid 
N bottom stage 


Superscript 


b boiling point 

Cc critical condition 

k iteration number 

fo) pure component, or ideal state 


R reduced state 
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APPENDIX 3-1 TEST PROBLEM NUMBER - 1 


*x*** PROBLEM STATEMENT **%*x 


TYPE OF COLUMN 1 
NO. OF PLATES 8 
NO. CF COMPONENTS 4 
PRESSURE AT TOP PLATE (PSIA) 300.00 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 2 
ENTERING PLATE OF FEED 1 1 
FEED QUANTITY 0.200008 02 
PRESSURE OF FEED 300.00 
TEMPERATURE OF FEED (F) 90.00 
ENTERING PLATE OF FEED 2 8 
FEED QUANTITY 0.100008 03 
PRESSURE OF FERED 300.00 


LIQUID RATIO, OF FEED 0.0 


*** RATE CF CONVERGENCE WITH 


6] SUM OF SGUARES OF RESIDUALS 


TEMPERATURE 
Oe2SS960E 03 0.55960E 
0255960E 03 0.55960E 
VAPOR FLOW RATE 
02.85000F 02 0.90000F 
O.90000E C2 0,90000E 


0.90000E 
0.95000E 


1 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 

0,2 58296E 03 0.59013E 
0.SS8020E 03 0.57014E 
VAPOR FLOW RATE 

0,2 S5961E 02 0.,91427E 
O2.92951E 02 0,.93739E 


z SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0.257689E 03 0.58159E 
O0.57837E 03 0.57080E 
VAPOR FLOW RATE 
0286080F O02 0.91825E 
0.293598E O02 0.941 &86E 


3 SUM OF SQUARES 
TEMPERATURE 
Owm579LI1E O03 0.S58466E 
O.S7750E O03 0.56956E 
VAPOR FLOW RATE 
O0.86012F 02 O0.91674E 
0.93388F O02 0.94026E 


03 0.58311E 
03 0.55555 


02 0.92603E 
C2 0.9S54C4E 


QF RESIDUALS 


G3 
03 


0.56538 1e 
0,55443E 


02 
02 


0.92284E 
0.95387E 


ITERATION NO, 


3 aE aK 


0.32670E-01 


03 
03 


0.55960E 


02 
02 


0.90000E 


0. 36086E-03 


03 0.59047E 
03 


O02 0, 92212E 
02 


0.14544EF-63 


03 0.58328E 
OS 


02 0,.92957E 
02 


0.76246FE-05 


03 
03 


0.58484E 


02 
02 


0.92 708E 


03 


02 


03 


02 


03 


02 


03 


C2 


10] 


0.55960E 


0.90000E 


0.58676E 


0.92497E 


0.58199E 


0.93245E 


0. 58216E 


0.93000E 
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03 
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RRR COMPUTED RESULTS **xxx 
NUMBER OF ITERATIONS = 
COLUMN PRESSURE (PSIA) O©.30000E 03 


PRESSUR 54 


TEMPERATUR » ENTHALFY AND LIQUID FRACTION OF FEED 
0.300F O3 PSIA O.900E 02 F 0.585EF 04.857... 1,000 
0. 300E 03 PSIA 0.558E F 0.432E 06,.8,T..U. 0.40 
eX STAGE VARIABLES *xx 
PLATE NO. TEMPERATURE LIG. Fe. Re VAP. F. Re HEAT DUTY 
F LB WOL/TIME LB MOL/STIME BTU/T IME 
1 O11 SSPreEresa OL 25662E 02 0. S60PZ2E OGL) COSC 
= O.12506E 03 0.26372E 02 O.91674E 02 0.0 
3 O0.12621E 03 0.26696E 02 .0.92384E 02 0.0 
4 0.12524E 03 02 26988E 02 OS BLTUSE (020.0_0 
2 GLTS256F 03 .0. 27376E 102 9.0 993000 2 102 1.0/6 
6 Oe-11790E C3 0. 280165E 02 0.93388E 02 0.0 
Fg 9~10996E 03 0.2°9374E 02 0.94028E 02 0.0 
8 0.94829E 02 0.339886 (02 ». 02952878 (02 /.620 
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*** MATERIAL BALANCES * xx 


FEEDS 
PLATE NO. 1 PLATE NO. 8 
LB MOL/STIME MOL FRC. LB MOL/TIME MOL FRC. 
ME THANE 0.0 0.0 0.70000E 02 0.70000F 00 
E THANE 0.0 c.0 0.15000F 02 0.15000F 00 
PROPANE 0.0 0.0 0.10000E 02 0.10000E 00 
NBUTANE 0.0 0.0 0.40000E O01 0.40000E-01 
NPENT ANE 0.0 0.0 O.10000F 01 06.10000F-01 
OCTANE O0.20000E 02 0.10000E 01 0.0 0.0 
TOTAL 0.20000E 02 O.10000F 03 
PRODUCTS 
TOP PRODUCT BOTTOM PRODUCT 
LB MOL/TIME MOL FRC. LB MOL/TIME MOL FRC. 
METHANE 0.67326E 02 C.78275E 00 0.26796F 01 0.7883SE-01 
E THANE O0.12514E 02 0.14549F 00 0.24845E 01 0.7309SE-01 
PROPANE 0.54286E 01 0.63114F-01 0.45680F 01 0.13440E 00 
NBUT ANE 0.30626E 00 0.35606F-02 0.36935 01 0.10867F 00 
NPENTANE 0,184156-03 0.21409E-05 0.99981F 00 0.29417E-01 
OCTANE 0.43764E 0CO 0,50881E-02 0,19562F 02 0.57557E 00 
TOTAL 0.86012E 02 0.33988E 02 


6o=-"0 O01 0 
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td Gqeave4d.o _ So 351064.6 
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APPENDIX 3-2 TEST PROBLEM NUMBER - 2 


¥¥%%x* PROBLEM STATEMENT ** xx 


TYPE, OF COLUHN 1 
NO. OF PLATES 6 
NO. OF COMPONENTS 15 
PRESSURE AT TOP PLATE (PSTA) 1310.00 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 2 
ENTERING PLATE OF FEED 1 1 
FEED QUANTITY 0.362508 02 
PRESSURE OF FEED 1310.00 
TEMPERATURE OF FEED (?F) ~§ 90 
ENTERING PLATE OF FEED 2 a) 
FEED QUANTITY 0.11809E 04 
PRESSURE OF FEED 1310.00 
TEMPERATURE OF FEED (F) 45.00 
NO. OF INTER COOLFRS OR HEATERS 1 
PLATE NC. OF INTER COOLER OR HEATER 1 6 


HEAT LOAD OF INTER COOLER OK HEATER -0.58000E 06 


= : 
7 
- 
OM.OFLt 
thst 
ones 
a f 


9 GUeONE.O- 


*%** RATE OF CONVERGENCE WITH ITERATION NO. 


0 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 
0.47960E 03 0.48860E 03 0.48460E 
O247960E 03 
VAPCR FLOW RATE 
0.93986E 03 0.95000FE C3 0.97500E 
O211000E 04 


1 SUM OF SQUARES CF RESTDUALS 
TEMPERATURE 
0.49032EF 03 0.48878E 03 0.48781E 
0.48296E 03 
VAPOR FLOW RATE 
0.92373E 03 0.96418E 03 0.98453E 
O-10811E 04 


= SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

0.48915EFE 03 0.48606E C03 
0248522E 03 

VAPOR FLOW RATE 

0.93488E 03 0.994S99E 03 
0211052E 04 


O0.4E416E 


0.10153E 


3 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 
02.49032E 03 O.48771EF 03 0.48538E 
0.248367E 03 
VAPOR FLOW RATE 
0.92874E 03 0j99307E C3 O0.J101I5I1E 
0.211023E 04 


4 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

O0.49077E 03 0.48985E 03 0.48828E 
02.48397E 03 
VAPOR FLOW RATE 

O0~93056E 03 0.99617F 03 0.10169E 
O02 10929E 04 


3 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 
0.48973E 03 0.48883E O23 
0O.48413E 03 
VAPOR FLOW RATE 
Oe. 93157E 03 0.99866E O03 
0210912E 04 


0.48762E 


O-1GIS7E 


ae a 
0.-89902E-01 


03 0.48460E 
03 0,10000EF 


0.85200£-02 


03 0.448625E 
O03 0.210045E 


0,25760E-02 


03 0.248322cE 
04 0210347E 


0.7768835-—03 


03 0.48370E 
04 0.210353E 


0 .89358E-04 


03 0.48€90E 
04 0,10351€ 


0.52666EF-04 


03 0,48672E 


04 0.10376E 


O03 


O4 


04 


03 


04 


03 


04 


04 


03 


04 


105 


0.48460E 


0,10500E 


0.48496E 


0210374E 


0,48367E 


0. 10644E 


0 .48304E 


0.10627E 


0. 105785 


0.48584E 


0.10585€E 


03 


04 


03 


04 


03 


04 


03 


04 


03 


04 


04 


La 7 


vu 


sO 


a1) 4008.0 


13006720 


befueyr.f) £0 


heave, a0 


sche.) £4 


7m ©, ii 2 
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VL awe EO 


ee 
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Ss oe 
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- ro. 2 


Welule #0 
s()= (SRE CR.O 
sauTanaam at 
iehhngt 6 vagaaene EO aev@ne.0 EO BFTOCAL0 
to 3YCLS*.0 
Sao wou? BodaV 
1400 t.take wala ye =O WIeee.e eo Sa208e,0 
#0 3052001 .0 
b= MNase 60 2 anuaEeS +o exnauen a mua 2 
SOUT ARI9M3T 
itvodjug 26 skayaned co 7086.0 £0 SETeRS.O 
Pasa Ween 
a0 boo no SERED Ae ep "0 3N21€%.0 
| ae] #0 3SIe0!.0 


6 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

0.49016E 03 0.48901E 03 0.48745E 
0.48384E 03 
VAPOR FLOW RATE 

0293008EF O03 0.99697E 03 0.10181E 
O210904E 04 


7 SUM OF SQUARES GF RESIDUALS 
TEMPERATURE 

O.49050E 03 0.48953E 03 0.48799E 
0248380E 03 
VAPOR FLOW RATE 

0.92995E 03 0.99648E 03 0,10172E 
0210895E 04 


0. 86600E-05 


03 0.448623F 


04 0.10361E 


0.50215£-05 


03 0.48662E 


04 0,10349EF 


03 


04 


04 
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0.48516E 


O~10571E 


0.10S557E 


03 


04 


a i ———— S * eee > > 2 


ie et 


dine ay Tae -_ 7s 


£0 DPSede.0 60 sdne.0 5, Me 


e%e*E* COMPUTED RESULTS ***xx 


NUMBER OF ITERATIONS 


COLUMN PRESSURE (PSIA) 


PRESSUR + 


O0.2131F 04 PSIA 


O,2131F 04 PSIA 


TEMPERATUR 


2 


ENTHALFY AND LIGUID FRACTION 


*** STAGE VARIABLES *** 


PLATE NO. 


On kwh 


TEMPERATURE 


F 
0.308S8E 
0.29931E 
0. 28386E 
0.27016E 
0.25649F 
0.24198E 


Ce 
02 
02 
C2 
02 
02 


LIG. Fs 


O2lZ1OOE O04 


Re 


LB MOLSTIME 


0. 1C278E 
0.12346E 
Osl41i6eE 
0.161S9E 
0.19582E 
0.28717E 


03 
03 
03 
03 
03 
03 


0.4505 02 F 0.342E OF. BST 550. 


VAP. F. 


Re 


LB MOL/TIME 


Oe J2Q995E 
0,.99648E 
0.210172E 
0, 10349E 
0. 10557 & 
0. 1OS9SE 


03 
03 
04 
04 
04 
04 
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OF FEED 
1,000 


0.0 


HEAT DUTY 
BTU/T IME 


0,0 
0.0 
0.0 
0.0 
0.0 
0.5 


aa 


Os | 
40 MOORE .O+ a0 2 


7 
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**e* MATERIAL BALANCES *x*x 


FEEDS 
PLATE NC, 1 PLATE NO, ral 
LB MOL/TIME MOL FRC, LB MOL/STIME MOL FRC. 
NITROGEN 0.0 0.0 0.20670F 02 0.17504E-01 
co2 0.0 Cc. 0 0.15250E 02 0.12914£-01 
METHANE 0.0 0.0 0.87215E ©3 0.73857E OO 
E THANE 0.0 0.0 0.17678E 03 0.14970CE 00 
PROPANE 0.0 0.0 0.68680EF C2 0.58161F-01 
IBUTANE 0.0 0.0 On.72800F O01 0,61650E-02 
NBUTANE 0.0 0.0 0.13670F 02 O,11576€F-01 
IPENT ANE 0.0 0.0 O-19700E O01 0.16682E-02 
NPENTANE 0.0 0.0 0.19500F O01 0,16512E-02 
HEXANE 0.0 0.0 O.13100FE O01 0.11094EF-02 
HYPTH-1 O.70000E-01 0.19310F-02 0.84000E 00 0.71134E-03 
HYP TH-2 0.29000F 00 O.80000F-02 0.25000F 00 0.21171F-03 
HYPTH-3 0.12990E 02 06.35835E CO 0.50000F-01 0.423425-04 
HY PTH-4 0.20430E 02 0.56359EF 00 O.10000E-01 0.84683F-05 
HYP TH-5 0.24700F O01 0.68138F-01 0,10000F-01 0.846823F-05 
TOTAL 0.36250E C2 0. 11809E 04 
PRODUCTS 
TOP PRODUCT BOTTOM PRODUCT 
LB MOL/TIME MOL FRC. 1B MOL/TIME MOL FRC. 
NITROGEN 0.19462E 02 0.20928E-01 0.12353F 01 0.43015EF-02 
co2 0.10284F C2 C.e110S9E-C1 0.49539F 01 6.17251EF-01 
METHANE 0.76128E O03 0.81862F 00 0.11101F 03 0.3865€F 00 
E THANE O0.11237E 03 6.12083E 00 0.64314F 02 0,22396F 00 
PROPANE 0.,.24443E 02 0.26284E+-01 0.4419646F 02 0.153SCE 00 
IBUTANE 0.10196E 01 0O.10964E-02 0.62578E 01 0.21791£-01 
NBUTANE 0.81845E 00 0.88010F-03 0.,12849F 02 0.44745£-01 
IPENTANE 0.14287E-01 0,1536€4E-04 O0.219557F 01 0.68101£-02 
NPENTANE 0.36481F-02 0,.39229F-05 0.19463F O01 O0.67776E-02 
HEXANE 0.25829F-04 0.27774E-C7 0213100F 01 0.4561 7E-02 
HYPTH-1 0.67593E-02 0.72684F-05 0.90320F 00 0.31452E-02 
HYP TH-2 0.10142F-01 0.109C6E-04 0.52979F 00 0.18449F-02 
HYPTH-3 0.15050EF 0O 0,16184EF-C3 0.412888E 02 0.448681E£-01 
HYP TH—4 0.90845E-01 0.97688E-04 0.20348F 02 0.70858E-01 
HYP TH-5S 0.46795E-02 0,.50320F-05 0.24753E 01 0.86195E-02 


TOTAL 


0.92995E 03 


0.28717E O02 
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APPENDIX 3-3 TEST PROBLEM NUMBER - 3 


**¥x** PROBLEM STATEMENT *%% xx 


TYPE OF COLUMN 2 
NO. OF PLATES 10 
NO. OF COMPONENTS 5 
PRESSURE AT TOP PLATE (PSIA) 300.00 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 2 
ENTERING PLATE OF FEED 1 1 
FEED QUANTITY 0.300008 02 
PRESSURE OF FEED 300.00 
TEMPERATURE OF FEED (Ff) 90.00 
ENTERING PLATE OF FEED 2 5 
FEED QUANTITY 0.10001E 03 
PRESSURE OF FEED 300.00 
LIQUID RATIO OF FEED 0.0 
SPECIFIED TOP PRODUCT 0.950008 02 
NO. OF LIQUID SIDE STREAMS 1 
LEAVING PLATE NO. OF STREAM 1 9 


QUANTITY OF SIDE STREAM 0.15000E 02 


a 


wl ee 


**%* RATE OF CONVERGENCE WITH 


0 SUM OF SQUARES GF RESIDUALS 
TEMPERATURE 

0.255960E 03 0.60960E 
0267960E 03 0.73960E 
VAPOR FLOW RATE 

Ow» 9SSO000E 02 0,.10500E 03 
0.75000E O02 0.50000E 


0.60S60E 
0.80960E 


O.10S00E 
0.S0000E 


1 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

0,~57872E 03 0,60439E 
O0.67867EFE 03 O0.74414E 
VAPOR FLOW RATE 
O.95000EFE O02 0.10354E 
O21 37500EFE 02 0.39664E 


03 
03 


0,60685E 
0.81110E 


03 
02 


0.1 OS2S5E 
024 8696E 


2 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

0.2458692E 03 0.60066E 
0.268078F O03 0.74897E 
VAPCR FLOW RATE 

0. 9SG00E O02 0,10S75E 
Q.18750E 02 0,.585052E 


03 
03 


0.60527E 
0.81 056E 


03 
Ce 


0.10603E 
O.48S61E 


3 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

O258748E 03 0,.60044E 
O.n69075EF 03 0.74993E 
VAPOR FLOW RATE 
O295000E 02 0.10404E 
O.15721E 02 0,38944E 


03 
03 


0. C0S36E 
0.80665E 


03 
02 


O0210€45E 
0.49323E 


4& SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

0.-58766E 03 0.60060E 03 0.60509F 
02.69482E 03 0.74976EFE C3 0.80EF73E 
VAPOR FLOW RATE 

06. S5000FE 02 010399F 03 0.10639F 
0.16389E 02 0.41181F 02 0.49545E 


5 SUM OF SGUARES OF RESIDUALS 
TEMPERATURE 

02e58767E 03 0.60069E 03 

02.69494E 03 0.74959E C3 
VAPOCR FLOW RATE 

O0.95000E 02 04~-10400E 03 
0.16141E 02 0,40574E 02 


0,60f544E 
0.80845E 


0.-10639E 
0.49229E 


ITERATION NO. 


3K IK Be 


03 0.6096€0E 
03 0.85960E 


03 0.11000E 
02 0.,50000EF 


0.19667E 00 


0, GO396E 
0.87 721E 


0.10830E 
0.51 238E 


0.44551E-01 


0.60 269E 
0.884496 


0.10796E 
D.SI277E 


0,122246-02 


03 
6 


0.603775 
0,.885605F 


03 
02 


0.210@13E 
0.50426E 


0.25923E-04 


03 0.60277E 
03 0.88565E 


03 0.10807E 
02 0.50775E 


st 51 2k -05 


03 
03 


0.60347E 
0.885505 


02£10805E€ 
0, 50 536E 


03 
02 


03 


03 
02 
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0.60960E 
0.953960E 


C211 S00E 
GC.65000EF 


0.59480E 
0.°94438E 


OelLIBie 
0.65246E 


0,.58921E 
0.94654E 


O..11026E 
0.64463E 


0.-S5S9213E 
0.294735E 


O.11014E 
0.64129E 


GO, o9085e 
0.94753E 


0.11022E 
0. 64046E 


0.59113 
0.94755E 


O211008E 
0.63 740E 


03 


03 


03 
02 
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*KKK COMPUTED RESULTS ***kx 


NUMBER OF ITERATIONS D 


COLUMN PRESSURE (PSIA) 


PRESSUR » TEMPERATUR 4,4 


02300E 03 PSIA 

0.300E 03 PSIA 

*** STAGE VARIABLES *** 
PLATE NC. TEMPERATURE 

F 

0 SL2ECTER OS 
Oe-14109E O83 
O.14584E C3 
0.14387E 02 
Oe TS PSS 0s 
0.23534E 03 
0.28999E 03 
0.348685E 03 
0.42590E 03 
0.-48795E 03 


ouMoUOMONO OH FW Nh & 


~~ 


ENTHALFY AND 


0.900EF 02 F O.877E 


O-446E O02 F 0.4166 


2 ee) ie 


0.30000E 03 


LIQUID FRACTION 


04 B.T. WV. 


06 -B.T 2. 


Re 


LB MOL/STIME 


0.38997E 
0.41392E 
0.43050E 
0,45077E 
0. 51TTSYE 
0. 7S584E 
0.84239E 
0. 8S546E 
0. BE2TSoe 
0. 20010EFE 


02 
02 
02 
02 
02 
02 
02 
02 
02 
02 


VAP, me 


Re 


LB MOL/S TIME 


0. 95000E 
0. 10400E 
0. 10639E 
0,.10805E 
O.11008E 
O.16141E 
O2.40574E 
0.49229E 
2 DU SSIOE 
0.O3740E 


02 
03 
03 
03 
03 
02 
02 
02 
02 


02 


111 


OF PED 


1.000 


0,0 


HEAT DUTY 
BTU/STIME 


89B4E O06 


7 a . 


eee 


re 


' - 


) @Gaevesia Ge 1nbVEe sO 


*%x* MATERIAL BALANCES 


FEDS 


METHANE 
E THANE 
PROPANE 
NBUTANE 
OCTANE 
TOTAL 


PRODUCTS 


ME THANE 
ETHANE 
PROPANE 
NBUTANE 
OCTANE 
TOTAL 


PLATE NOG, 


6 MOL/STIME 


e 
0.0 
0.0 
0.0 
0.0 
0Oe-30000E O02 
0.30000E 02 


° 


BK OK 3c 


i 
MOL FRC, 
0.0 
0,0 
0.0 
0.0 
O.100COE O1 


TOP PRODUCT 


LB MOL/STIME 
0.79998E 02 
O,66644E 01 
0.62250E O01 
O.15469E O1 
0.56545E 00 
0.9S5000E O02 


LIQUID SIDE STREAMS 


METHANE 
ETHANE 
PROPANE 
NBUTANE 
OCTANE 
TOTAL 


MOL FRC. 
0.84209E 00 
0. 7O1S1E-—01 
O.6552%o-—01 
0, 162836-—01 
Cp 9DS21b=— OZ 


LIQUID’ 4) 9} 


LB MCL/STIME 
0.1586 7E-02 
0.47338 (2-02 
0.31347E C0 
0.32002E O1 
0.11480E 02 
0.15000E O02 


MOL FRC, 
0. 10576E—03 
6, 3159 1-03 
C.208S8E-01 
0.21334E 00 
O0.76534E 00 


We 


PLATE NO. S 


iB MOLSTIME 
0.80000E 02 
0.66700E O01 
0.66700E O1 
0.66700E 01 
0.0 

0. 10001F 063 


BOTTCM 


LB MOL/STIME 
0,31 999E-03 
0. 894 08E- 03 
GO. AS31435 60 
0219230E 01 
0,17S9S54E 02 
O0.20010E 02 


MOL FRC. 
0.79992E 00 
0.66693E-01 
0.66693E-01 
0.66693E-01 
0.0 


PRODUCT 

MOL FRC. 
0.21599 2E-04 
0.44682E-04 
026564 4E-02 
O0.9610 15-01 
0.83727E 00 


a 04 STALL 


ows 4 


OO ASCOT RO BOO0OS.0, 
16-36 0086.0 


10-70 7O00,.0 
19-2 OA4656 
0.0 
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APPENDIX 3-4 TEST PROBLEM NUMBER - 4 


*eKRK PROBLEM STATEMENT **k&x 


TYPE OF COLUMN 2 
NO. OF PLATES 16 
NO. OF COMPONENTS 15 
PRESSURE AT TOP PLATE (PSTIA) 21D OO 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 2 
ENTERING PLATE OF FEED 1 1 
FEED QUANTITY 0.26590EF 03 
PRESSURE OF FEED 270.00 
TEMPERATURE OF FEED (F) “6.00 
ENTERING PLATE OF FEED 2 6 
FEED QUANTITY 0.17108E 04 
PRESSURE OF FEED 270.00 
TEMPERATURE |OF FEED (F) 50.00 
SPECIFIED TOP PRODUCT 0.67700E, 03 
NO. OF INTER COOLERS OR HEATERS 1 
PLATE NO. OF INTER COOLER OR HEATER 1 1z 


HEAT LOAD OF INTER COOLER OR HEATER O.40000E O07 


*** RATE OF CONVERGENCE WITH 


0 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0253450E 
0.4 59460E 
0. 65960E 
0..70960E 


03 
03 
03 
03 


0.S54S60E 
0.60960E 
0.€6960E 


VAPGR FLOW RATE 


0.67 700E 
Ow SS000E 
0. SSOO00E 
0210000E 


03 
03 
03 
04 


0. 75000E 
0 .90000E 
Os. 1LOO0OE 


C3 
vs 
03 


03 
Oz 
04 


0.5S5S60E 
0.61960E 
0.67S60E 


0.80000EF 
0.90000E 
0.10000E 


1 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0.2 5S3689E 
O.53914E 
0,.60433E 
O.-70095E 


03 
03 
03 
03 


0O.S54821E 
0.55369E 
06627 TSE 


VAPOR FLOW RATE 


0267700EF 
Ow D1L497E 
04 1Z2337E 
0.291 728E 


03 
03 
04 
03 


0.86131E 
0,.84562E 
0. 130335 


c3 
03 
03 


03 
GS 
C4 


0.S5S2¢7e 
0 .S651LTE 
0.64326E 


0. 8Y279E 
0.935518E 
6. 75189E 


= SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


Oe SIS30E 
02 S3709E 
0.57983E 
0.69805E 


03 
03 
03 
03 


0 »HD4344E 
0.54796E 
02.61055E 


VAPOR FLOW RATE 


0.2 67700E 
Ow 92805E 
0. 93463E 
O.84437E 


03 
03 
03 
03 


0. 27268E 
0. 86192E 
0, 97750E 


03 
03 
03 


03 
03 
o3 


0.54865E 
0.55249E 
0. 62ES7E 


0.2€870CE 
0294227E 
0.5 785715 


3 SUM OF SQUARES GF RESIDUALS 
TEMPERATURE . 


025338 8E 
Oe S3785E 
02 57763E 
02 69326E 


03 
03 
03 
o3 


0,.54432E 
0.54904E 
0.60168E 


VAPOR FLOW RATE 


On O67 700E 
0. 9513 8E 
Ow LOBG6E 
02 76309E 


03 
03 
04 
03 


0 .88606E 
0.90112E 
Oe LIQ20E 


03 
03 
03 


03 
03 
C4 


0.9O0f44E 
0.10153E 
0,67154E 


ITERATIGN NO. 


ok 


0.873582" 06 


03 
03 
03 


03 
03 
04 


0.56 S60E 
0.63460E 
0.268960E 


0.8S5000E 
0. 90000E 
O210000E 


0.25738S5E-01 


03 
03 
03 


03 
03 
03 


0.55564E 
0.257 869E 
0.63804E 


0. 86833E 
0.10020E 
0.11 S63E 


0.19695E-01 


03 
03 
03 


3 
03 
03 


0. 5526S5E 
0255755 
0. G29S7E 


0, 88636F 
0.94431E 
0, 84886E 


0245025EF-02 


03 
03 
03 


03 
04 
03 


0.554236 
0. 55ST1E 
0. €1904E 


0.2.91234E 
0.10506E 
0.65781€ 


0S 
03 
03 


03 
03 
04 


03 
03 
03 


03 
03 
03 


03 
03 
o3 


03 
O04 
03 
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0.57960E 
0.64960EF 
0.69960E 


0.90000E 
0.900C00E 
O0.10000E 


0.554684E 
0.59154€E 
0.04958E 


O.87146E 
0.11198E 
0.94992E 


« SS293E 
0.56533E 
0 -64423E 


0. 89035E 
0.93353E 
0.83915€ 


eSS417E 
0.56619E 
0.63512E 


0.91963E 
0.10681E 
0.76298E 


03 
O3 
03 


03 
03 
04 


03 
03 
03 


03 
04 
03 


03 
03 
03 


03 
03 
03 
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€o 
Ee 
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bi 


a0 
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ea 


0 
Pt 


co 


y oy 
ego 
2%) 


to 
o 
tf 
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300008 .0 


2000628,0 
390000 .0 
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geuee?.o 
o42@107,0 
AO OD. 0) 


Fee tw, 
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TAGeee, 6 


WVU? .O 
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94,0000 ,0 
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ca 
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42 EAR yt 


Te fa €p 
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e0 
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ED S683 19.0 


by 


fO-S80401.0 a9 2eRaloe 2O wUe 
wet 6ea.0 Oh 
see te2.6. 20 
ee ta Ww Bb 


7 


wr-inw &o 
Cane, CO 
iV4e3A,9 £60 


S0—-4 2EO01te, 0 
W502 O-K8 


THVH St CO 
Ue ee 2) 


4a SUM GF SQUARES OF RESIDUALS 


TEMPERATURE 
0.2.53 727E£' O03 
0.-S3779E O03 
0.S57565E 03 
0.69376E 03 


0.54612E 
0.54892E 
0.60058E 


VAPOR FLOW RATE 


Oa 67700E 
0. SS5S55E£ 
0.2 LLI6O0E 
GO. 80236E 


03 
03 
04 
03 


0.8&8162E 
0.90842E 
O.-11409E 


03 
03 
04 


0. €9S34E 
0.1025€E 
0.7168S5E 


=a SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0.53636E 03 
S.53783E 03 
0.57643E 03 
0.69369E 03 


0.54700E 
0.54895E 
0.60185E 


VAPOR FLOW RATE 


Oe267TOOE 
0.95486E 
OO. LL 105E 
0.4 79319E 


03 
03 
04 
03 


0. 88356E 
0-9O0746E 
0.2 11302E 


03 
03 
03 


C3 
03 
04 


0.55226E 
0 5S3S7E 
0.60S37E 


0.90209E 
0.210252E 
0.7C228E 


0.267 03E-03 


03 
03 
03 


03 
04 
03 


0.55534E 
0.S5775E 
0.61 849E 


0.90761E 
0.10637E 
O.77481E 


0.213340E-04 


03 
03 
03 


03 
04 
03 


0.555 13E 
0.55790E 
0. 61933E 


0.90991E 
0210639E 
0.75446E 


03 
03 
03 


03 
04 
03 


03 
OS 
03 


03 
04 
03 


1S 


0. SS369E 
0.56433E 
0.263432E 


O.91S72E 
0.10857E 
0.83737E 


2 S5404E 
0,.S56476E 
0263468E 


0.91958E 
0,10848E 
O.31767E 


03 
03 
03 


03 
O4 
03 
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XR COMPUTED RESULTS **%%% 


NUMBER OF ITERATIONS 


COLUMN PRESSURE 


PRESSUR , 


TEMPERATUR 


(PSIA) 


= 


ENTHALFY AND LIGUID FRACTION 


O.27CCOE 03 


0.270EF 03 PSIA -0,500E ©1 F-0.915E 06 B.T..U.~ 


Owz270E 03 PSIA 


*** STAGE VARIABLES *x* 


PLATE NO, 


of Gi fo 


ON O 


TEMPERATURE 


fF 
0, 76756E 
0.87 396E 
0.92660E 
6.95531E 
0 .94436E 
0.78234E 
0. 89355E 
0.93972E 
0. 98299E 
0.210516E 
0.11683E 
0214225E 
O.14977E 
05 1S973E 
0.17508E 
0.23409E 


02 
02 
02 
C2 
02 
02 
C2 
02 
C2 


ciG. Fe 


Re 


(6B WOL/S TIME 


0.47246E 
0.490S9E 
0,.4S881E 
0, 508485 
0. S4376E 
0.22072E 
O. 23249E 
0 .23636E 
0. 23845E 
0.24102E 
0.24299E 
0. 20C2Z0E 
0.20542E 
OS 2LIZ4e 
0.209295 
0212997E 


03 
03 
03 
03 
03 
O04 
04 
O04 
04 
O04 
04 
04 
04 
O4 
04 
04 


O.-SOO0E 02 F-0.146E 06.8.7..U. 


VAP. F. 


Re 


LB MOLSTIME 


0.67700E 
0, 88356E 
0.90209E 
0. GOSSIE 
0,91958E 
0.95486E 
0. 90746E 
0. 102525 
0. 10639E 
0. 10848E 
OSLLTOSE 
0. 11 302E 
0. 70228E 
0. f5446E 
0. S1V67E 
02/9319E 


03 
03 
03 
C3 
03 
03 
03 
04 
04 
04 
04 
04 
03 
03 
03 
03 
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1.000 
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eee MATERIAL BALANCES **x 


FEEDS 
PLATE NO, 1 PLATE NO. ra 
LB MOL/STIME MOL FRC, LB MOL/TIME MOL FRC. 
NITROGEN 0.0 0.0 0.700C0E 00 0.4091 7E-03 
co2 0.0 0.0 0.12400F 02 0.72481E-02 
ME THANE 0.0 0.0 0.16740E 03 0.9784S9F-01 
ETHANE 0.0 0.0 O.47440E 03 0.2773CE 00 
PROPANE 0.0 0.0 0.44020F 03 0.25731E 00 
I BUTANE 0.0 0.0 O.64300E 02 0.37585F-01 
NBUTANE 0.0 0.0 0.12880E 03 0,.75286F-01 
IPENTANE 0.0 0.0 0.19500E 02 0.113986-01 
NPENTANE 0.0 0.0 0.19400F 02 0.11340E-01 
HE X ANE 0.0 0.0 0.13100F 02 0.76572E-02 
HYP TH-1 O.SO0000E 00 0.18804E-02 0.90000F 01 0,52607E-02 
HYPTH-2 O.21000E 01 O0.78S77E-C2 0.530C0E 01 0.30980E-02 
HYP TH-3 0.95300E 02 0.35841E 00 0.12830F 03 0.74994E-01 
HYP TH-4 0.14990£ 03 C.56375E 00 0.20330F 03 0.11883E 00 
HYP TH-5 O.18100£ 02 0.68071E-01 0.24700EF 02 0.14436F-01 
TOTAL 0.26590E 03 0.17108E 04 
PRODUCTS 
TOP PRODUCT BOTTOM PRODUCT 
LB MOL/TIME MOL FRC. LB MOL/TIME MOL FRC. 
NITROGEN 0.69924E 00 (C.10328E-02 0.62390F-11 6.48003E-14 
co2 0.12402E 02 0.18319F-01 0.39272E-02 0.30216€E-05 
METHANE 0.16734F 03 0.24718E 00 0,12452E-04 0.95808EF-08 
ETHANE 0.44615E 03 0.65901F 00 0.28300F 62 0.21774F-01 
PROPANE 0.49821F 02 0.73591E-01 0.39038E 03 0.3003€F 00 
I BUTANE O.28764E 00 C.424E7F-03 0.64012E 02 0.49252E-01 
NBUTANE 0.12321E 00 6.18200F-03 0,12868F 03 0.99005F-01 
IPENTANE O0.17947E-03 0.26509E-06 0.19500F 02 0.150036-01 
NPENT ANE 0.42297E-04 0.62477E-C7 0.194C0E 02 0.1492€EF-01 
HEXANE 0.52087E-07 0.76928E-10 0.13100£ 02 0.10079E-01 
HYPTH-1 0.71291E-02 C.10520E-04 0.94929F 01 0.7303SE-02 
HYP TH-2 0.84500E-02 0.12481E-04 0.73916€ 01 0.56871E-02 
HYP TH-3 0.10585E 00 0,.15635E-03 0.22349F 03 0.17196E 00 
HYPTH-4 0.52168E-01 C.7705GE-04 0.35315F 03 0.27171£ 00 
HYP TH-5 0.23296E-02 0.34410F-05 0.42798E 02 0.32929F-01 
TOTAL 0.67700E 03 0.12997E 04 
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APPENDIX 3-5 TEST PROBLEM NUMBER - 5 


**x** PROBLEM STATEMENT **%%% 


TYPE OF COLUMN 2 
NO. OF PLATES 20 
NO. OF COMPONENTS 10 
PRESSURE AT TOP PLATE (PSIA) 464.70 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 2 
ENTERING PLATE OF FEED 1 1 
FEED QUANTITY 0.118308 04 
PRESSURE OF FEED 464.70 
TEMPERATURE OF FEED (FP) 5.00 
ENTERING PLATE OF FEED 2 17 
FEED QUANTITY 0.12792E 04 
PRESSURE OF FEED 464.70 
TEMPERATURE OF FEED (F) 50.00 
SPECIFIED TOP PRODUCT 0.642208 03 
NO. OF INTER COOLERS OR HEATERS 2 


PLATE NO. OF INTER COOLER OR HEATER 1 
HEAL LOAD OF INTER COOLER OR HEATER 


PLATE NO. OF INTER COOLER OR HEATER 2 
HEAT LOAD OF INTER COOLER OR HEATER 


m 
-O0.84000E 06 


7 
ae. 109508 OF 


118 


Bt 


OY . 08 


v 
at &O0QUK,0- 


y 
i POR VOT oe 


***x RATE OF CONVERGENCE WITH 


0 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


02 45960E 
O.247460E 
0.49960E 
0. 53960E 


03 
03 
03 
O3 


0. 46460E 
0.46960E 
0.50960E 
0. 54G960E 


VAPOR FLOW RATE 


0-2 64220E 
0.2 92000E 
0.12000E 
02.13500E 


03 
03 
04 
04 


0265000E 
0. 93000E 
0.212500€E 
Oe14000E 


0.46S60E 
0.47160E 
0.2.51960E 
0.57S60E 


0265000E 
0.1 1OO0E 
O.1l25S00E 
0.25000E 


1 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 

O02 45565E 03 0,45264E 
On.47942E O03 O47719E 

0.49002E 03 0.49070E 

0.48891E 03 0.48960E 
VAPCR FLOW RATE 


0. 64220E 
0.214205E 
0218985E 
0. 16981E 


re SUM OF SQUARES OF 


03 
04 
04 
O4 


TEMPERATURE 


O.47113E 
0247395E 
02.48557E 
0249 740E 


03 
03 
03 
03 


0.75424E 
0, 182715 
0.18769E 
0. 16766E 


0.4648S5E 
0.46816E 
0.48726E 
0.49973E 


VAPOR FLOW RATE 


04 64220E 
0.212630E 
O2lVO9I7TE 
0.215 796E 


03 
04 
04 
04 


0.74350E 
O213500E 
0.16955E 
0.15570E 


03 
03 
03 
03 


03 
04 
04 
04 


R 


03 
03 
03 
03 


03 
04 
04 
O4 


0.~45348E 
0.4&€38E 
0.48944E 
0.52538 


0.86727E 
0.1 7827E 
0.18228E 
0.39196E 


ESTIDUALS 


0.45828E 
O.47827E 
0.48926E 
0.S52S577E 


0.81€44E 
0.159 78€ 
O2.16702E 
0.24016E 


a SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 
O42.47916E 03 0,48656E 
Oe47783E 03 0.47154E 
O249009FE 03 0.49145E 
Oe.49688E 03 0.499046E 
VAPOR FLOW RATE 


Ow 64220E 
GO. LVISSE 
Ow 1S546E 
0. 15269E 


03 
04 
04 
04 


G-« 74532E 
0.12195E 
0.15540E 
06152325 


03 
63 
03 
Os 


03 
04 
04 
C4 


O.47E74E 
0.48206E 
0.4S283E 
0.SZ2ESGE 


0,80132E 
0.14358E 
0.15466E 
Os.Z2 OS9TE 


ITERATION NO, 


03 
03 


03 
04 
04 
03 


es 


0.46460£F 
0.4796 0E 
0.52960E 
0,58 960E 


O0.70000E 
0.11500E 
0.13000E 
0.35000E 


0.552325 00 


03 
O03 
03 
03 


03 
04 
O04 
03 


0,45153E 
0.49046EF 
O.48978E 
0.54 740E 


0210128E 
0.18770E 
0217860E 
0.41545E 


0.10724E 060 


03 
O03 
03 
03 


03 
O04 
04 
03 


0.244 73SE 
0.4821 8E 
0249202E 
0. S4373E 


0.91110E 
0. 16887E 
0.16419E 
0.439416 


0.29842E-01 


03 
03 
03 
03 


03 
O04 
04 
03 


0.45667E 
02.48642E 
0.-49431E 
0.54137E 


0.85894E 
0.15194E 
0.15406€E 
0.246 084E 


03 


03 


vg 


0.46960E 
0,48960E 
0.53460E 
0,60960E 


0.90000E 
O212000F 
0.13000E 
0 s400G6E 


0,47 323E 
0249124E 
0.4888S8E 
O,S8497E 


O.12901E 
O.19170E 
O« PTSYSE 
0 .43686E 


0.46689E 
0.48430EF 
0249472E 
0.58428E 


O0.21124€E 
0.17158E 
0.16076E 
Os49770E 


0.46999E 
0.48865E 
0.49565E 
0.58321E 


O.10132E 
02 15496E 
O0.e.15318E 
0.52409E 


03 
03 
03 
03 


03 
O04 
04 
os 


03 
03 
03 
03 


04 
04 
04 
03 


03 
03 
03 
03 


04 
04 
04 
03 


03 
03 
03 
03 


O04 
04 
04 
03 


V0OC 96.0 
206086,.0 
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304004,0 
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4 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0.47500EFE 03 
Oe247803E 03 
O249000E 03 
0.49784E 03 


0.48528E 
0.47023E 
O0,49161E 
O,49971E 


VAPOR FLOW RATE 


0.64220E 
0.210752E 
0214882E 
0.214887E 


03 
04 
04 
04 


O.75101E 
Os 1 PSE&G6GE 
0.14913E 
0. 14900E 


03 
03 
03 
03 


03 
04 
04 
04 


0.48S938E 
O.4E&108E 
0 .49316E 
0.52706E 


0.81421E 
0. 1 3623E 
O.14899E 
O.18134E 


+) SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
02.47584E 03 
0.2.47903E 03 
0.48924E 03 
Ow 49765E 03 


0.48274E 
0246920E 
0.49088E 
0.49958E 


VAPOR FLOW RATE 


0.64220E 
0.104986 
Os. 14503E 
0. 14688E 


03 
04 
04 
04 


0.75026E 
O.1TP2iSE 
O214562E 
0.14726E 


03 
03 
03 
C3 


03 
04 
04 
04 


O.48569E 
0.,.47S98E 
0.49251E 
0.52621E 


0.81698E 
0213161E 
0214591E£ 
0.16946E 


6 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0247562E 03 
0248088E 03 
O248918E 03 
O0.49749E 03 


0248354E 
0 .46964E 
O.49077E 
0,4993SE 


VAPOR FLOW RATE 


0.264220E 
O02 10418E 
0. 14267E 
02.14567E 


03 
O4 
04 
O04 


0. 75098E 
0.2.11036E 
0214346E 
0.214625E 


C3 
03 
C3 
03 


03 
04 
04 
C4 


O0.4E€557E 
0.47986E 
0.49239E 
0.52563E 


0. 81634E 
O.12876E 
0.14398E 
0.16€203E 


rs SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
O247571E 03 
0.48185E 03 
0.48930E 03 
O0.49744E 03 


0.48342E 
0O2.47044E 
0.49088E 
0.49932E 


VAPOR FLOW RATE 


O0264220E 03 
0210388E 04 
02141266 04 
02414496E 04 


0.75120E 
0.10962E 
O. 14213E 
0.14567E 


0.4ESS2E 
0.48016E 
0.49246F 
0.52527E 


0.817335E 
Ge t27Z22F 
0.14285E 
0.15 754E 


0.8971 76-02 


a 
03 
03 
03 


03 
04 
O04 
03 


0.46866E 
0. 4857 6E 
0.49480E 
0.54070E 


0.86066E 
0, 14472E 
0.142894E 
0.46384E 


0.24243E-02 


O3 
03 
O23 
03 


03 
04 
04 
Os 


O,47314E 
0.48487E 
0.49424E 
0, 54 C29E 


0.87 022E 
0.1402S9E 
0214619EF 
0,459 75E 


0. 72298E-035 


03 
03 
03 
03 


03 
04 
04 
03 


0.47223E 
0.48477E 
0.49410E 
0. 54000EF 


0,86911E 
0.13746E 
0.14449F 
0. 45405E 


0.21769E-035 


03 
O03 
03 
03 


03 
04 
O04 
03 


0.47244E 
0.4849 1E 
0.49412E 
0. 53981E 


0. 86943E 
0.213578E 
0.14 349E 
0.4491SE 


C3 


120 


0 .47363E 
0.248828E 
0.49636E 
0.58294E 


0.97460E 
O0.14790E 
C.14879E 
6.53012E 


O.4778BE 
O.48747E 
0.49597E 
0,58270E 


0.97044E 
0.14373E 
O,14647E 
0.530416 


0e47972E 
0.48739E 
02.49583E 
0.58252E 


O.97055E 
0.14112E 
0.14503E 
0.52890€E 


0.47S94E 
0248752E 
0.49582E 
0.2.58240F 


0. 96886E 
0. I SOS6E 
Ow~144186E 
0.527425 


03 


03 
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can 


8 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0.47570E 
0.-48215E 


0.48938E 03 


Os 49745E 


03 
O03 


03 


0.248352E 
0.47086E 
0O.490G4E 
0.49S31E 


VAPOR FLOW RATE 


0.2.64220E 
0. 10372E 
Oe 14045EF 
0. 144555 


03 
04 
04 
O4 


0. fSUSTE 
0.10930E 
O.1l4146E 
0. 14533E 


9 SUM OF SQUARES OF R 
TEMPERATURE 
O0.47568E 03 0.48347E 
0248233E 03 0.47103E 
O248940E 03 0.49096E 
On 49745EFE 03 0,49931E 
VAPOR FLOW RATE 


02 64220E 
0.10369E 
Oe 14000E 
0.214432E 


10 


03 
04 
04 
04 


TEMPERATURE 
0~47568E 03 0,48346E 
O02 48244E 03 0,47113E 
0.48942E 03 0.490S97E 
02e49745F 03 0.49930E 

VAPOR FLOW RATE 


02 64220E 
O21lD371E 
02 139746E 
O.14419E 


03 
04 
04 
04 


O.75147E 
O.10917E 
0-14105E 
0. 14513E 


0. 75152E 
0.-10913E 
Oe 14082E 
0.214502E 


03 
03 
03 
03 


03 
04 
04 
C4 


03 
O3 
03 


03 


03 
04 
04 
04 


O.4E5S9E 
0.48042E 
0.49251E 
0.S2SC08E 


0.81762E 
O.12E47E 
0.14220E 
O0.15495E 


ESIDUALS 


0.4 E600E 
Oe48054E 
0.49253E 
02524S9E 


0,81739E 
0212610€ 
O.14164E 
0.153828 


SUM OF SQUARES OF RESIDUALS 


O.-4E597E 
0.48060E 
02.49254E 
0-52494E 


0,&1799E 
0.12592E 
O.14164E 
0,15273E 


0.65867E-04 


03 
03 
03 
03 


03 
04 
04 
03 


On4T727GSE 
0.48504E 
0.49416E 
0.53 970E 


O~87014E 
0,13486E 
0214292E 
0.44565E 


0221433E- 04 


03 
03 
03 
03 


03 
04 
O04 
03 


0.472S90E 
0248512E 
Oe49417E 
0. S3964E 


0, 87053E 
02-13437E 
0.,14260E 
0. 4432S9E 


0.66766E-05 


03 
03 
03 
03 


03 
O04 
04 
G3 


0.4728S9E 
0.248515E 
0.4941 7E 
0.539606 


0, B7074E 
0. 13412E 
0.14242E 
0.44 183€ 


03 
03 
03 
03 


03 
C4 
04 
Oe! 


03 
03 
03 
03 


03 
04 
04 
03 


03 
03 
03 
03 


03 
O04 
O04 
C3 
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0.48020E 
0.48700E 
0.49584E 
0. SB233E 


0.96879E 
0.13866E 
0.14369E 
0.526275 


0.480415 
0.,48764E 
0.49585E 
0.58229E 


0.,96930E 
O0.13817E 
0.14 342E 
0.52549E 


0.48048E 
0248766E 
O.49584E 
0.58227E 


0,.96965€E 
0.13790E 
0.14326E 
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**EEE COMPUTED RESULTS *#%%x 


NUMBER OF ITERATIONS 


COLUMN PRESSURE (PSIA) 


PRESSUR 


0.2465FE O03 PSIA 


0246S5E 03 


TEMPERATUR 


TA 


10 


ENTHALFY AND LIQUID FRACTION 


**%* STAGE VARIABLES **x 


PLATE NO, 


OANA ESE WD & 


10 


TEMPERATURE 


F 
0,16084E 
0.238635 
0.26366E 
0.13289E 
0.20876E 
0222844E 
oO, LI SSiTe 
0,21005E 
0, 25S5352E 
0. 28063E 
0.29817E 
0.31365E 
0.,32936E 
Oe34574E 
0, 36244E 
0.37852E 
0,39698E 
0.65341E 
0.80002E 
0.12267E 


02 
C2 
02 
C2 
G2 
02 
C2 
02 
02 
C2 
02 
C2 
C2 
02 
02 
02 
02 
02 


03 


Des ae 


0.46470E 03 


C.SQ0E C2 F 0.485E O07 


Re 


LB MOL/STIME 


0. 12923E 
0. 1S5@8E 
0214115E 
0. 15104E 
0. 15779E 
0. 16321E 
0. 1E000E 
0, 18820E 
0.19198E 
0. 1S382E 
0.19490E 
0. 1SS7e2E 
0. 1S6S50E 
0219734E 
0. 19827E 
0.1S9910E 
0219728E 
0. 226135 
0.232450E 
0. 18200E 


04 
04 
04 
04 
04 
O04 
O04 
o4 
O04 
04 
O04 
04 
O4 
04 
04 
04 
04 
04 
04 
04 


O.500E 01°F-0.949F 07.847. 50. 


B.1..2U. 


VAP, F. 


Rs 


LB MOL/ TIME 


0.64220E 
GO. 21 see 
0,.81799E 
O.8S7074E 
0. S696S5E 
O.10371E 
0. i10913E 
0. 12592E 
0.213412E 
QO, 1STS0E 
0. 13974E 
0.14082E 
0, 14164E 
Oe 14242E 
0. 14326E 
0,144195 
Q0.214502E 


0.15273E 


0.441836 
0.52497E 


O3 
03 
03 
03 
03 
O4 
04 
04 
O04 
04 
O4 
G4 
o4 
04 
04 
04 
O04 
03 
03 
03 
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OF FEED 


1,000 


0.0 


HEAT DUTY 
BYIUSTIME 
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0.0 
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—-0.84000E 06 
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—-Q.10950E O07 
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*e* MATERIAL BALANCES **x 


FEEDS 


METHANE 
E THANE 
ETHYLENE 
PROPANE 
HYDROGEN 
BENZENE 
TOLUE NE 
O-XYLENE 
PROPYLEN 
ECHEXANE 
TOTAL 


PRODUCTS 


METHANE 
ETHANE 
ETHYLENE 
PROPANE 
HYDROGEN 
BENZENE 
TOL VE NE 
O—-XYLENE 
PROPYLEN 
ECHEXANE 
TOT AL 


PLATE NO. 


LB MOL/STIME 


O9BOE 02 
S660E 02 
S150E o2 
3660E 03 
1830E 04 


» 
+ 
» 
» 


0.0 
oO, 0 
0.0 
0.0 
0.0 
0.f73261VE 03 
0.7 
0.2 
0.5 
Owe 
0.1 


1 


MOL FRC, 


G 
0 
je 
0.0 
0.0 
C.67CCOE 00 
0.60000E-01 
C.20000E-01 
0.500C0E-O01 
0.20000EF 00 


TCP PRODUCT 


LB MOL/STIME 
0.11822E 03 
0.81437E oc 
0.54739E 02 
0.16102E-11 
0.45792E 03 
0.66477E 00 
0.16202E-01 
0.14592E£-02 
0.97940E 01 
0.28427E-01 
0.64220E 03 


MOL FRC. 
C.18409E 00 
©. 12681E-02 
UsoO esrTe- oF 
0.25074E-14 
0. 71 30SE 00 
0.10385146-02 
0. 25229E— 04 
0. 227% 226-05 
W.152516-o1 
C.44266E-04 
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PLATE NO. 17 


LB MOL/STIME 
0.12792E 03 
0.27887E 03 
0O.4131SE 03 
Os.12790E 01 
0.45796E 03 


2792E 04 


BOTTCM 


LB MOL/TIME 
O.97011E O1 
C.22TCQSE oS 
0, 35845E 03 
O.12790E 01 
0. 37136E- 01 
O.f91S4E 03 
O.70964E O02 
GO. 236538E 02 
0.49 356E 02 
0.23657E£ 903 
0218200F 04 


MOL FRC, 
0.9999SE-01 
O0.21800EF 00 
0.32300EF 00 
0.9998 2E—-03 
0,.35800E 00 


PRODUCT 

MOL FRC. 

0.53302E-02 
0.15278E 00 
0.196S5E 00 
0.70274E-03 
0.20404E-04 
0.43512E 00 
0.38991£-01 
0.1299SE-01 
0.2711SE-01 
0.12998E 00 
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APPENDIX 3-6 TEST PROBLEM NUMBER - 6 


xe PROBLEM STATEMENT ***%% 
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TYPE OF COLUEN 2 
NO. OF PLATES 30 
NO. OF COMPONENTS lig 
PRESSURE AT TOP PLATE (PSIA) 475.00 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 3 
ENTERING PLATE OF FEED 1 1 
FEED QUANTITY 0.18200E 04 
PRESSURE OF FEED 475.00 
TEMPERATURE OF FEED (F) -4.00 
ENTERING PLATE OF FEED 2 15 
FEED QUANTITY 0.28900E O04 
PRESSURE OF FEED 475.00 
TEMPERATURE, OF FEED (F) -4.00 
ENTERING PLATE OF FEED 3 22 
FEED QUANTITY 0.29600EF 04 
PRESSURE OF FEED 475.00 
TEMPERATURE OF FEED (F) ~4.00 
SPECIFIED TOP PRODUCT 0.15200E 04 


eee RATE OF CONVERGENCE WITH 


0 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0.45960E 
0.49460E 
0.52960E 
0. 56460E 
02SS960E 
0» 63460E 


03 
03 
03 
03 
03 
03 


0.46660E 
0.50160E 
0.53660E 
0.57160E 
0.60660E 
0.64160E 


VAPOR FLOW RATE 


O02 1S200E 
0. 1E000E 
02. 20000E 
0. 20000E 
04 24000E 
O-16000E 


04 
04 
04 
04 
04 
04 


0.16000E 
0.18S500E 
0.20000E 
0.20000E 
0.2S5000E 
0.18000E 


03 
03 
03 
03 
03 
03 


04 
C4 
04 
04 
C4 
04 


0.4 7360E 
0~.SCEBCE 
0 .54360E 
0.57860E 
0.61360E 
0.64860E 


O.216S5C00E 
O0.19C000E 
O.20000E 
0.21000E 
0.10C000E 
0.20000E 


1 SUM GF SQUARES OF RESIDUALS 
TEMPERATURE 


0.42869E 
0. 45075E 
02.47645E 
0. 54059E 
0. 57584E 
0.64831E 


03 
03 
03 
o3 
03 
03 


0.41278E 
0.45959E 
0.4834SE 
0.55060E 
0.58386E 
0.64682E 


VAPOR FLOW RATE 


O02 1S200E 
0217944E 
0.20838E 
0.23568E 
O02 29701E 
0. 22480E 


04 
04 
04 
04 
04 
04 


0.15998E 
O0,18565E 
0.21004E 
0224323E 
0.31196E 
0.271 73E 


03 
03 


0.41360E 
0.46543E 
0.4°S496E 
GO, 55853 
0.60938E 
0.653 24E 


0. 1€455E 
0.19236E 
0.21048E 
0.25783E 
0.18933E 
0.25940E 


a SUM GF SQUARES OF RESIDUALS 
TEMPERATURE 


0.46041E 
Oe248790E 
026 48262E 
0253222E 
0.5604 8E 
0.61968E 


03 
03 
03 
03 
03 
03 


O02«45808E 
0.48265E 
0.48751E 
0254416E 
0.56357E 
0e62819E 


VAPOR FLOW RATE 


0. 15200E 
Ow 18649E 
0.221092E 
0.27690E 
0. 36507E 
0.2.29314E 


04 
04 
04 
04 
04 
04 


0.16475E 
0.19286E 
0.21219E 
0. 29746E 
02«38315E 
0.27681E 


03 
O3 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 


0.46S4SE 
024 7S34E 
0.49475E 
OP apes War fz 
0 .S59861E 
0,62936E 


0.16742E 
0. 1SE64E 
0.21336E 
02.31806E 
02.2€44S9E 
0.32042E 


ITERATION NO, ¥*x 


O-12434E O1 


03 
03 
OS 
03 
O3 
03 


04 
04 
04 
O4 
04 
04 


0.48060E 
0. 515605 
0.55060E 
9... 55 5605 
0. €20E60E 
0.65560E 


0.17000E 
O219SO00E 
0.20000E 
0222000E 
0.12000E 
0.22000E 


0.63660E CO 


03 
03 
03 
03 
03 
03 


O04 
04 
04 
04 
O4 
04 


O4 


0242466E 
0.46°38E 
0.S50993E 
0.56400E 
0. 64057EF 
0.65715E 


0. 16899E 
0, 19942E 
O.20947E 
o.27011€ 
0.1469 1E 
0.27930E 


0.48466EE 
0.478S6E 
0.50383E 
0.5555 7E 
0.60959E 
0,€4257E 


O.17211E 
0.2041 7E 
0.21469E 
Oe.33491E 
0.22408E 
0.30160E 


03 
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0.48760E 
0. 52260E 
0,55760E 
0.59260£ 
0.2E2760E 
0.,66260E 


O0.17500E 
0.20000E 
0.20000E 
0.23000E 
O0.14000E 
0.24000E 


0.24385S9E 
0.47295E 
0.52506E 
0,56951€E 
0.63660E 
02.f6747E 


0.,17385E 
0, 20662E 
0.20813E 
0.28365E 
0.2.23436E 
0.27 892E 


O.-49091E 
0.48004E 
O.SVITSE 

» 558i5€E 
O.61971€E 
02.662 09E 


O.179IGDE 
0.2.20943E 
0.21768E 
0.35014E 
0.23490E 
Oe3SIVIE 


03 
03 
03 
03 
OS 
03 


04 
04 
04 
04 
04 
04 


03 
03 
O3 
03 
03 
03 


04 
O04 
04 
04 
04 
O04 


03 
03 
03 
03 
03 
03 


O04 
04 
04 
04 
04 
04 


s0at64,.0 &O 
so4gne.0 
s05v68.0 (© 


£0 


go04908@.0 
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3 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0249 160E 
0» 49046EF 
0.249826E 
0251577E 
0. 54409E 
0. 58385E 


03 
03 
03 
03 
03 
03 


0 249409E 
0 -49474E 
0.S50062E 
O.53041E 
0.53595E 
0.-S8447E 


VAPOR FLOW RATE 


Oe 1S200E 
0.20592E 
0. 2083 3E 
Ow 32006E 
0244219EF 
0. 36496E 


04 
04 
04 
04 
04 
04 


0. 18081E 
0 4.20824E 
O220874E 
0.2.36789E 
0245542E 
0.40788E 


C4 
04 
O04 
04 
04 
04 


0.18536E 
0.20826E 
0.210S1E 
0,.3SE&87E 
0.32675E 
Os3F601E 


4 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0.47125E 
02 50609E 
0.251982E 
0.251782E 
0.54830E 
0. SE0944E 


0.15200E 
0. 19686E 
0.20489E 
0. 28814E 
0. 40895E 
0242425E 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 


0.50790E 
0.51048E 
0.52084E 
0.53080E 
0.53910E 
0. 5S9499E 


VAPOR FLOW RATE 


0.19129E 
0,20207E 
0.20619E 
0.32941E 
0.4193SE 
0.43260E 


03 
OS 


04 
04 
O04 
04 
C4 
04 


0.52800E 
0.51575E 
0.252147E 
0,.54C78E 
0.56759E 
0.60459E 


0.19533& 
0.20435E 
0.20798E 
0.36163E 
0.26286E 
0.43729E 


5 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


02 46463E 
0, S2i157E 
0252712E 
0. 5234 8E 
0.2 54880E 
0.2.58949E 


03 
03 
03 
03 
03 
03 


0.47893E 
0.51864E 
02.52779E 
0, 53278E 
0.53992E 
0.59507E 


VAPOR FLOW RATE 


Ow 1 S5200E 
0-2 20096E 
Oe 2O841E 
Oe 2923 1E 
02.40789E 
Ow 41 797E 


04 
04 
04 
04 
04 
04 


0. 18838E 
0.20303E 
0. 20935E 
0.232759E 
0.41871E 
02.42676E 


03 
03 
03 
03 
c3 
03 


04 
04 
04 
04 
04 
04 


0.50751E 
60.52055E 
0.S2801E 
0.54124E 
0.56789E 
0.,60477E 


0.19S62E 
0.20559E 
0.21 087E 
6.35888E 
0.25889E 
0243136E 


0.31891E-01 


oS 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 


O- 


OS 
03 
03 
03 
03 
03 


04 
04 
04 
04 
O04 
04 


0. 47234E 
O.49471E 
0. 50429E 
0254423E 
0.57410E 


0. 61 SS5£ 


0, 19058E 
0.2081 0EF 
0.221462E 
0.4151 2E 
0.36599E 
0.425146 


15235E-61 


» S2283E 
0.5184 0E 
0.2.52010€ 
0. 54677E 
0.57 892E 
0.62345E 


0219270E 
0, 2044 SE 
0.21 086E 
0. 38398E 
0,37323E 
C.43917E 


0.86465E-02 


03 
03 
03 
03 
03 
03 


04 
O04 
04 
04 
04 
04 


GaSZOLTZE 
0, 523768E 
0. S2592E 
0.54695E 
0.S57S09E 
0.62367E 


0. 20413E 
0.20719EF 
0.21341E 
0.38198 
0.36539F 
0.43435E 


03 


03 
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0.47994E 
0.49591E 
0.49725E 
0.54580E 
O.S7446E 
0.65613E 


0,19887E 
0.20811E 
O.22411E 
0.243205E 
0.379 09E 
0.41 766E 


O.5S1034E 
O0.51925E 
0,5083S5E 
0.54939E 
0. 58521E 
0,65867E 


0.19263€ 
0.20438E 
0.21793E 
0.39884E 
0.40634E 
0243992 


0.52911€E 
0.2.52603E 
0.51276E 
0.54968E 
0.58S515€E 
0.65886E 


0.201 89€ 
0.20789E 
O.21991E 
0. 329742E 
0.40058E 
Oe43771F 


03 
03 
03 
OS 
03 
03 


04 
04 
04 
04 
04 
O04 


03 
03 
03 
03 
03 
03 


O04 
04 
04 
04 
04 
04 


O03 
03 
03 
O03 
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03 


O04 
04 
04 
04 
O4 
O04 
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6 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0. 46315E 
0.52201E 
0. S2646E 
O.S2710E 
0.54917E 
0. 58948E 


03 
03 
03 
03 
O03 
03 


O.46411E 
0.51924E 
0.52852E 
0.S3569E 
0.54036E 
0.59506E 


VAPOR FLOW RATE 


0.215200E 
0.21088E 
0. 21386E 
0.29693E 
Os 407416 
O241587E 


04 
04 
04 
04 
04 
O04 


0.187S58E 
0.20975E 
0.21441E 
0.33024E 
0.41818E 
0.242467E 


O.47751E 
O.S1854E 
0.52974E 
0.54271E 
0.56813 
0.60480E 


0.19926E 
O0.21C69E 
0.21535E 
0.35885E 
0.25724E 
0.42S56E 


F i SUM GF SQUARES OF RESIDUALS 
TEMPERATURE 


0.46898E 
0.50231£ 
Oe.51741E 
0. 53027E 
0.54920E 
0258912E 


O3 
03 
03 
03 
03 
03 


0.46634E 
0.508496 
9 .52207E& 
0.53848E 
0.5406CE 
0.59468E 


VAPOR FLOW RATE 


0.15200E 
0.2.21919E 
0. 22146E 
0... 30556E 
0.40943E 
0.241503E 


04 
04 
04 
04 
04 
04 


0.18566E 
0.21974E 
0.2.22201E 
0.33812E 
0.41993E 
0.42363E 


03 
03 
C3 
03 
03 
03 


C4 
04 
04 
C4 
04 
C4 


0.46674E 
O.SO0S31E 
0.52€4CE 
0 .54428E 
0.5681 7E 
0260442E 


0.19530E 
0.21912E 
0.222248E 
0O.3E4C1E 
0.25815E 
0.42E864E 


8 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0447824E 
0.248395E 
0.504455 
Ow S3203E 
0254909E 
0. 58848E 


03 
03 
03 
03 
03 
03 


0.48413E 
0.49324E 
0.50803E 
0.S4106E 
0.S4067E 
0.59397E 


VAPOR FLOW RATE 


O-. 1S200E 
0.20526E 
0.22281E 
0. 3228 SE 
0.241343E 
Ow 41 431E 


04 
04 
04 
04 
04 
04 


O.18407E 
Oe21411E 
0.22529E 
0.35547E 
0.242326E 
0.42266E 


C3 
03 
03 
03 
03 
03 


04 
04 
04 
04 
C4 
04 


0e.4€202E 
0.S50079E 
0.51403E 
0.54633E 
0.56811E 
0.60369E 


0.18844E 
0.21E72E 
0.22778E 
0.3774SE 
0.26C083E 
0242768E 


0.62034E-—02 


O3 
03 
03 
03 
O3 
03 


04 
O04 
04 
04 
04 
04 


0.5008S5E 
0.520€1£ 
0.52843£ 
0.54765E 
0. 57920E 
0.62372E 


0.20820E 
O.Z212°07E 
02.21727E 
0.38127E 
0. 36234E 
0.43305E 


0 -43366E-02 


03 
03 
03 
8S 
03 
03 


04 
04 
O04 
04 
O4 
04 


Oe 


03 
03 
03 
03 
03 
03 


04 
O4 
04 
04 
04 
04 


O.47441F 
0.51 026E 
0. 52804E 
0.54824E 
O.S7904E 
0.62335E 


0.20453E 
0.21955E 
6.22368E 
0.38450E 
0.36155E 
0.2.43235E 


17647E-02 


0.4788SE 
0.50322E 
0.52002E 
0. S4923E 
0.57869E 
0. 62263E 


O.19147E 
0.21987E 
0.23054E 
0.39325E 
0. 36157E 
0.43152E 


03 
03 
03 
03 
03 
03 


04 
04 
04 
O4 
04 
04 


o3 
03 
03 
O3 
03 
03 


04 
O04 
04 
04 
04 
04 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 
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0.518S7E 
0.S52371€E 
0.51570E 
0.55008E 
0.58512E 
0.65888E 


0.21234E 
@.Z1325E 
0.22301E 
0,39677E 
0.39870E 
0.243750E 


0.488575 
G.SPSi0£ 
0.51809E 
9. 558020E 
0.58478E 
0.65855E 


G.2U3 72E 
0.22057E 
0.2284°E 
6. 3991SE 
0.39823E 
Os.43751E 


0247886E 
0. 50339€ 
0.51632E 
0.55040E 
0.S5S8419E 
0,65798E 


0.219678E 
0.22088E 
0.23638E 
0.40475E 
0.39815E 
Ow43714E 


03 
03 
03 
03 
03 
03 


O4 
04 
04 
04 
04 
04 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 


O3 
03 
O03 
03 
03 
03 


o4 
04 
O4 
04 
o4 
04 
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9 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
O.47671E 03 
02.48972E 03 
0. 50765E 03 
O0.52978E 03 
0.S4933E 03 O.5S4074E 
OeS8S851E 03 O.59401E 

VAPOR FLOW RATE 
O0215200E 04 0.18542E 
O.20179E 04 0.20813£ 
Oe21867E 04 0.22029E 
O.32159E 04 0.3567SE 
Ow41359E 04 0.42303E 
O241457E C4 0.42301E 


0.48743E 
02.49406E 
0.5098e&E 
0.540 17E 


10 
TEMPERATURE 
O.47476EF 03 
Oe49840E 03 
0.251084E 03 
0.252763E 03 
0.54946E 03 0.54081E 
& S8S861E 03 0.5941 1E 
VAPOR FLOW RATE 
O.-15200E 04 0.18612E 
O.-20101E£ O04 0.2050S5E 
O221615EFE 04 0.21740E 
On 31784E 04 0.35537E 
O-41377E 04 0.42293E 
O0241439E 04 0.42286E 


0.48545SE 
0 .49899E 
0. 51 286E 
0.53875E 


ti 
TEMPERATURE 
02.47330E O03 
0.-50504E 03 
0251302E 03 
0. S2653E 03 
0454949E 03 0.54085E 
02e58869E 03 0.59420E 
VAPOR FLOW RATE 
O.15200E 04 0. 18631E 
Oe 20327E 04 0.20554E 
02.21526F 04 0.21654E 
0.31388E 04 0.35224E 
On 41 367E 04 0-42279E 
0.41421E 04 0.422735E 


0.48194E 
0.50491E 
G.SIS4tE 
0.537f1E 


03 
03 
C3 
03 
o3 
03 


c4 
04 
04 
c4 
04 
04 


O04 
04 
04 


SUM OF SQUARES OF R 


03 
03 
03 
os 
C3 
03 


C4 
04 
04 
04 
04 
c4 


0.2492C08E 
0.4S988E 

S136SE 
0.54628E 
0.5€818E 
0.60375E 


0.1S066E 
6.21381 
0.222S55E 
0.37SSS5E 
0.2608S5E 
0.42814E 


SUM OF SQUARES GF RESIDUALS 


0.49423E 
0.50187E 
0.515405 
0.54562E 
0.56826E 
0.60385E 


0.19282E 
0.20924E 
0.21S926E 
@.2T97SE 
0.26C67E 
0.42602E 


ESIDUALS 


0.49068E 
0.505615 
0.51769E 
0.54497E 
0. 56832E 
¢..60395E 


6.1939¢E 
0 220832E 
O.21SL1e 
$.37 8656 
0.26047E 
0.42788E 


0.99069E-03 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
O04 
04 


0.49145E 
0,50416E 
0.51 760E 
0, 549575 
0.57874E 
0.62270E 


0.19296E 
0.21599E 
0.22593E 
0. 39504E 
G,. 361'S55E 
0.2.43205E 


0. S8286E-—s 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 


0.49905E 
0. 505S7E 
0.51 760E 
0.5S4S941E 
0.57 383E 
0.62282E 


0. 19576E 
0. 21 266E 
GO. 22252E 
0.39577E 
0. 36128E 
0.43197E 


0.32767E-03 


03 
03 
63 
03 
03 
03 


04 
04 
O04 
04 
04 
04 


Oe49874E 
0.50 767E 
0.51894E 
0.54912E 
0.57890E 
0.62291E 


0. 19830E 
0.21123E 
0.22095E 
0.3950 8E 
0.36101E 
0.4318SF 
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0O.48931E 
0.50629E 
0.51 306E 
0.S5077E 
0.58423E 
0.65806E 


0.19630E 
0.21 746E 
9.23319E 
0.40566E 

a39827E 
0.43772E 


0.49947E 

« 50865E 
O.SILI40E 
0.55065E 
0.S58432E 
0.65816E 


G.21TST9TE 
O.21483E 
0.23012E 
0.240622E 
0.39803E 
0.243772E 


0. 50370E 
0.51 038E 
0.51127E 
0.55078E 
0.58440E 
0.65824E 


0.20116E 
O.21361E 
0)..22820€ 
0.40600E 
0.39780E 
0.43 766E 
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Ha 


Sgescs 
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i2 
TEMPERATURE 
0.247290E O03 
0. 50586E 03 
0. 51398E 03 
0.52640E 03 
0.54943E 03 0.54085E 
0.58872E 03 0.59423E 

VAPOR FLOW RATE 
Oe1S200EFE 04 0.18620E 
0.20638E 04 0.20821E 
O221555E 04 0.4216S5E 
0O.231153E 04 0.34957E 
0241354E 04 0.442276E 
On41411E 04 0.442265E 


O0.47974E 
0.50789E 
Q0.51670E 
0.53722E 


is 
TEMPERATURE 
O0e47336E O03 
08 50237E 03 
Oe-SI1405E 03 
02 52684E 03 
0.54933FE 03 0.54080E 
0258870E 03 90.59420E 
VAPOR FLOW RATE 

O.15200E 04 O0J18S601E 
0O-20792E 04 0.21058E 
O0.21649EF 04 0.21786E 
Ow S1L1OZE 04 0.34825E 
O02. 41345E 04 0.42282E 
0e41409FE 04 0.422636 


0.47969E 
0.50659E 
a ..516ait& 
0.537225 


14 
TEMPERATURE 
0247392E 03 
0249946E 03 
9251382E OS 
Oe252740E 03 
02454925EFE 03 0.54075E 
0.58866E 03 0.59417E 
VAPCGR FLOW RATE 
O21 S200FE 04 0.18590E 
0.20752E 04 0.211 01E 
0e.21739E 04 0.421865E 
0.31170FE 04 0. 34814E 
0.41333E 04 0.42284E 
0041411E 04 0.42264E 


0.-48070E 
0,50390E 
0.51648E 
0.53 746E 


C3 
03 
03 
03 
C3 
03 


C4 
04 
04 
04 
04 
04 


03 
03 
03 
03 
03 
03 


C4 
04 
04 
C4 
04 
04 


03 
03 
03 
cS 
03 
03 


04 
04 
04 
C4 
04 
04 


SUM OF SQUARES OF RESIDUALS 


O0.4°641E 
0.50€&62E 
02£519628E 
0.54452E 
0.56833£ 
0.60398E 


0219411E 
0.20S89E 
0.21843E 
Os3761I4GE 
0.26€043E 
0.4277G9E 


SUM OF SQUARES OF RESIQUALS 


0.48484E 
0O.SO89TE 
0,.51S76E 
0.54431Ee 
0,.56831E 
0.6039SE 


Os427 785 


SUM OF SQUARES OF RESIDUALS 


0.48S561E 
0.SO738E 
0.51956E 
0.54434E 
0.56827E 
0.6039i1E£ 


6.19333E 
0.2 1340E 
O.22015E 
0.37438E 
0.26061E 
0.42775E 


0.17100E-—03 


03 
03 
03 
03 
03 
03 


o4 
04 
04 
04 
04 
04 


0.49396E 
0. 50964E 
0252040E 
0.548 82E 
0. 57892E 
0.62294E 


0.19937E 
0.21184E 
GC, 22093E 
0.239404E 
0,.36090E 
0.43183E 


0.900°1E-04 


03 
O63 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 


0.49043E 
9,51 OS9E 
0.52124E 
0.54860E 
0.57890E 
0. 62292E 


0.19908E 
0.213616 
Os 222 73E 
05393 22E 
0.3609 1E 
O.43177E 


0.46719E-04 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
O04 


0.49008E 
0.50982E 
0.521395 
0.548593 
0.57886E 
0. 622675 


RL 1TVSSTE 
0.21497E 
0.22252E 
0,39275E 
0.36095E 
0.43176E 


03 
O3 
03 
C3 
03 
03 


04 
O04 
04 
O04 
04 
04 
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0,S50107E 
0.51150 
0.51199E 
0.55063E 
0.S8443E 
0.65826E 


0.20354E 
0421 384E 
0,22767E 
0,405S56E 
0.39767E 
0243761E 


0.49659E 
0.51192€E 
Os 51 2E6E 
0.5S046E 
0.58440E 
0.65823E 


0.20395E 
0.621503E£ 
0.228i10E 
0.40518E 
0.39768E 
0.43 756E 


0249468E 

oD1L1L74E 
0.Si337E 
0. 55037£ 
0.58437E 
OO, Sasi: 


0.203 14€E 
0.21620E 
0.22880E 
0.40489E 
0.39769E 
0.2.43752E 
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1S 

TEMPERATURE 
0.47437E 93 
02. 49820E 03 
Ce oLseSe 03 
0252808E 03 
02 54919EF O03 0.54069E 
0.258863E 03 0.59413EF 
VAPGR FLOW RATE 
0.215200EF 04 0.18594E 
0420551E 04 0.20946E 
0.221830E 04 0.421941E 
O.231351E 04 0.34918E 
0241 305E 04 0.42270E 
0.41412E 04 0.442266E 


0 .248235E 
0.50108E 
0,515S4E 
0.53803E 


16 
TEMPERATURE 
0.47404E 03 
0.50021E 03 
0.251306E 03 
02 5279SE 03 
0.54925EF 03 O0.54071E 
0.58864E 03 + SY9415E 
VAPOR FLOW RATE 
0215200E 04 0.18604E 
O-. 20553E 04 0,20880E 
Oe.21784E 04 0.21913E 
02 31367E 04 0,34970E 
02441299F 04 0.422S59E 
0.41412E 04 0.42266E 


0.48198E 
0.50247E 
9. SLS9SE 
0.53810E 


03 
03 
03 
03 
03 
03 


C4 
04 
04 
C4 
04 
04 


03 
03 
03 
C3 
03 
C3 


04 
04 
04 
04 
C4 
04 


SUM OF SQUARES OF RESIDUALS 


0.488&20E 
0.50440E 
0.S51€95E 
0.54463E 
0.56822E 
0,60387E 


0.1S314E 
0.212856 
0.22086E 
0.37448E 
0.26059F 
O.42776E 


SUM OF SQUARES OF RESIDUALS 


0.48835E 
0.50484E 
0.51894E 
0,54476E 
0.S6824E 
0.60389E 


0619342E 
Osa twse 
0.22061E 
0.37495E 
O.26€C47E 
O0.42777E 


0.245701E-04 


03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
O04 


0.49253E 
O,50777E 
O.52101E 
0.54861£ 
0457883E 
0.62283E 


0. 197GiE 
0. 21542E 
0.223345 
0. 39241E 
0.35095E 
0.43176E 


0.23063E-04 


03 
03 
03 
03 
03 
03 


04 
04 
04 
O04 
04 
04 


0.49361E 
« 50 7S0E 
0.52 084E 
» S4873E 
0.57884E 
0.62285E 


0.19804E 
0.214336 
0.22308E 
0.39 265E 
0.36091E 
O.43177F 
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0.49560E 
0.5107 LE 
Os51L3S58E 
0.55034E 
0 .58434E 
0. 65817E 


0.20164E 
0.21712E 
0.22976E 
0,40450E 
0.39770E 
0.43 744E 


0Oe4A9748E 
0. 5)029E 
G.5E3S0E 
0.55043E 
0.58435E 
0.658198 


0.20198E 
8. 206326 
0.22958E 
0.40455E 
0.39770E 
0.43751€E 


03 
03 
03 
03 
03 
03 


04 
04 
04 
O04 
O04 
04 


03 
03 
03 
03 
03 
03 


04 
04 
04 
O04 
04 
04 
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RK COMPUTED RESULTS xx xx 
NUMBER OF ITERATICGNS 16 


COLUMN PRESSURE (PSIA) O.47500CE 03 


PRESSUR » TEMPERATUR », ENTHALFY AND LIQUID FRACTION OF FEED 
O.475E 03 PSIA -0.400E 01 F-C.618E 07 B.T..U. 1.000 
O2475E 03 PSTA -0.400E 01 F-0,805E O7 BiTeno. 1,000 
G247S5SE 03 PSIA -0.400E 01 F O04413E OF BLTSCU. 0,386 


#** STAGE VARIABLES **% 


PLATE NO. TEMPERATURE LIC. F. R. VAP. F. Re HEAT DUTY 
F LB MOL/TIME LB MOL/TIME BTU/TIME 
1 0.14441E 02 0.21604F 04 0.15200E 04 0.0 
2 0.22375E C2 0.22342E 04 0.18604E 04 0.0 
3 0.28752E C2 0.22804E 04 0.19342E 04 0.0 
4 0.34013F 02 0.23198EF 04 0.19804E 04 0.0 
cS) 0.37883E C2 0.23553E 04 0.20198F 04 0.0 
6 O.40609E C2 0.23880E 04 0.20553E 04 0.0 
7 0.42869E 02 0.24178E 04 0.20880EF 04 0.0 
6 0.45244E C2 0.24433E 04 0.21178E 04 0.0 
) 0.47902E 02 0.24632E 04 0.21433F 04 0.0 
10 0.S0686E 02 0.24784E 04 0.21632E 04 0.0 
11 0,53464E C2 0.24913F 04 0.21784E 04 0.0 
12 0.56333E 02 0.25061F 04 0.21913F 04 0.0 
13 0.59342E 02 0.25308E 04 0.22061F 04 0.0 
14 0.61243E C2 0.25S958E 04 0.22308F 04 0.0 
15 O.53705E 02 0.463267E 04 0.22958E 04 0.0 
16 0.68352E 02 0.66870E 04 0.31367E 04 0.0 
17 0.78503E 02 0.69395EF 04 0.34970F 04 0.0 
18 0.85165E 02 0.71165E 04 0.37495EF 04 0.0 
19 0.89130E 02 0.72355E 04 0.39265E 04 0.0 
20 0.908Z20E 02 0.473199E 04 0.40455E 04 0.0 
21 0.89645E 02 0.74159E 04 0.41299F 04 0.0 
22 0.81114E C2 0.&7547£ 04 0.42259F 04 0.0 
23 0.10864E 03 0.97591E 04 0.26047E 04 0.0 
24 O.11924E 03 0.10127E 05 0.36091E 04 0.0 
25 0.12475E 03 0.10291E 05 0.39770EF 04 0.0 
26 0.12904E 03 0.10377E 05 0.41412E 04 0.0 
27 0.13455E C3 0.104286 05 0.42266E 04 0.0 
28 0.14429F 03 0.10468F 05 0.42777E 04 0.0 
29 0.16325E 03 0.10525E 05 0.43177E 04 90.0 
30 0.19859E 03 0.61500F 04 0.43751F 04 0.30393E 08 
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*%%* MATERIAL BALANCES **x 


FEEDS 


HYDROGEN 
NITROGEN 
METHANE 
ETHYLENE 
ETHANE 
PRGPYLEN 
PROPANE 
I—BUTENE 
Crs26UTE 
NBUTANE 
NPENT ANE 
HEPTANE 
TOTAL 


HY DROGEN 
NITROGEN 
METHANE 
ETHYLENE 
E THANE 
PROPYLEN 
PROPANE 
I—BUTENE 
GIS2BUTE 
NBUTANE 
NPENTANE 
HEPTANE 
TOTAL 


PRODUCTS 


HYDROGEN 
NITROGEN 
ME THANE 
ETHYLENE 
FE THANE 
PROPYLEN 
PROPANE 
I-BUTENE 
CIS2BUTE 
NBUTANE 
NPENT ANE 
HEPTANE 
TOTAL 


PLATE NO. 


B MOL/STIME 


0.23660E C2 
0.19656E 03 
0.44408E 023 
O.77714E 03 
OesSIiZae “OS 
0.47320E 02 
O,18200E 04 


PLATE NC. 


LB MCL/STIME 
0.56240E 03 
0e20720E 03 
0.651205 505 
Os41440E 03 
O.41440E 03 
0.38480E C3 
O.11840E 03 
0.88800E C2 
0-14800E C2 
CO. S9200E O02 
6. 35520E a2 
O.88800E O1 
0.29600E 04 


0 
c 
6) 
O. 
C 
0 
GO, 1S 000E- 01 
C.10ECOE CO 
0.24400E 900 
C.427C0OE 00 


C.18200E 00 
C.26000E-G61 


22 

MOL FRC. 

C.2190G0E 00 
C.70000E-01 
G.22000E 00 
Ce~140CG0E 00 
0.14000E 00 
C~13C0CE 00 
C.40000E-O01 
©.300C0E-—01 
0.500C0E—-02 
0.20000E-01 
C.120C0E-01 
2 J0000E-O2 


TOP PRODUCT 


LB MOL/TIME 
0.56242E 03 
0.20720E 03 
0.67942E 03 
0.26266E 02 
0.29935E 00 
0,10018£-06 
0.26912E 01 
0.67539E O01 
0.99860E O01 
0.22224E 02 
0.26981E 01 
0.37131F-01 
0.15200E 04 


MOL FRC. 
C.37002E 00 
0.13621E 00 
0.446S99E 00 
0.172806-01 
0.19694E-03 
0, 65S05E-10 
0.17705E-02 
0.44423E-02 
0. 656¢S7E-02 
0,146€215£-01 
0.17 750E-02 
C.24428E-04 


PLATE 
iB MOL/STINE 
0.0 
0.6 


0.289CO0E 02 
0.57800E 02 
0214450F 03 
0.49130E O03 
O.31790E 03 
0.52020E 63 
0,17340E G3 
0.98260E 03 
02.15895E 03 
0.14450E O02 
0.2e@SO00E 04 


PEAT. 
iB MOL/STIME 


BOT TOM 


LB MOL/TIME 
0.3465 9E- 05 
0.66211E-03 
0.67292E 00 
0.44592E 03 
0.55860E 03 
0.87610F O03 
0.45727E 03 
0.79880F 03 
0.62229E 03 
O.17967E 04 
0.52301£ 03 
0.70613E 02 
0261500E 04 
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NO, iS 
WOL FRC. 
0.0 
0.0 


C.10000E-01 
0.20000E-01 
0.SO00C0E-01 
0.1700CE 00 
O.11000E 900 
C.1i800CE 00 
0.60000L-01 
0.240C0E 00 
0.5500CE-01 
0.5000G0E-02 


NO. 


MOL FRC. 


PRODUCT 

MOL FRC. 

0.56356E-09 
0.,1076€E-06 
0.1094 2E-03 
0.7250€E-01 
0.9082SE-01 
0.14246E 00 
0.74352E-01 
021298SE 00 
0.1011S€ 00 
0.29215E 00 
0.85042E-01 
0.11482F-01 
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APPENDIX 3-7 TEST PROBLEM NUMBER - 7 


wt PROBLEM STATEMENT ***+x 


TYPE OF COLUMN 3 
NO. OF PLATES 21 
NO. OF COMPONENTS 11 
PRESSURE AT TOP PLATE (PSTIA) 264.70 
PRESSURE DROP PER PLATE (PSIA) 0.0 
NO. OF FEEDS 2 

ENTERING PLATE OF FEED 1 7 

FEED QUANTITY 0.34500E 02 

PRESSURE OF FEED 264.70 

LIQUID RATIO OF FEED 0.0 

ENTERING PLATE OF FEED 2 13 

FEED QUANTITY 0.65500E 02 

PRESSURE OF FEED 264.70 

LIQUID RATIO OF FEED 1.0000 
REFLUX RATIO %A22 
LIQUID RATIO OF TOP PRODUCT 0.0 
SPECIFIED TOP PRODUCT 0.23000E 02 
NO. OF LIQUID SIDE STREAMS 1 

LEAVING PLATE NO. OF STREAM 1 4 

QUANTITY OF SIDE STREAM 0.15000E 02 
NO. OF VAPOR SIDE STREAMS 1 

LEAVING PLATE NO. OF STREAM 1 16 


QUANTITY OF SIDE STREAM 0.25000E 02 


0 SUM OF SQUARES GF RESIDUALS 
TEMPERATURE 


0. S3960E 
0.263210E 
Oe 72460E 
0.81710E 
0.290960E 


03 
03 
03 
03 
03 


0. 553 0CE 
O0.65060E 
0.74310E 
0. 2356 CE 


VAPOR FLOW RATE 


0.23000EF 
0. 94806E 
0.60306E 
0. 60306E 
0.85306E 


02 
02 
02 
02 
02 


0.S4806E 
0.94806E 
0,60306E 
0 .85306E 


0.S7660E 
0.669170E 
0.76160E 
0.8541 0E 


0.94806E 
0.60306E 
0.60306E 
0.85306E 


1 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0.53145E 
0.260341E 
0.268382E 
0. 76206E 
0.8656 2E 


03 
03 
03 
03 
03 


0.55484E 
0,62044E 
0.,70149E 
O.77724E 


VAPOR FLOW RATE 


0223000E 
0.87516E 
Oe 52847E 
0.44461E 
0.268851E 


G2 
02 
C2 
02 
02 


0.94806E 
0.85239E 
0. 51023E 
0 <GSS5S7£ 


03 
C3 
03 
03 


C2 
2 


02 
02 


0.56383E 
0.64183E 
0.7 29359E 
0.79410E 


0.9€213E 
0.52582E 
0.47090E 
0.6€701E 


2 SUM OF SQUARES OF RESIDUALS 


0.54503E 
02 S9B00E 
0266974E 
0. 72898E 
0.82006E 


VAPGR FLOW 


O02 23000F 
0.2 89762E 
0. 49858E 
0.44024E 
0.260559E 


TEMPERATURE 


03 
03 
03 
03 
03 


0. S6592E 
0.60635E 
0.68419E 
0. f3S501TE 


RATE 


02 
02 
02 
02 
02 


0, S4806E 
0 2-86482E 
02.48202E 
0.68264E 


0.57518E 
0. 63060E 
0.71615€ 
0.74434E 


0,.95819E 
0.50087E 
0,42324E 
0.667555 


*%** RATE OF CONVERGENCE WITH ITERATION NO. 


IKK 


0.87316E 00 


03 
oS 
03 
OS 


02 
02 
02 
02 


0.5951 0EF 
0.68 7/760E 
0.78010E 
0. &726 0E 


0,.9480€E 
0.60 306E 
0.60 306E 
0. 85306E 


0, 19655£° 00 


03 
03 
03 
03 


02 
02 
02 
02 


0.57266E 
0». €5746E 
0.73890E 
0.812945 


0. 94 526E 
0.53330E 
O.45771£ 
0.69375E 


0.19760E-02 


03 
03 
03 
03 


02 
02 
02 
02 


0.58235F 
0. 646405 
0.72049E 
0.75874E 


0.94950E 
0. 5021 2E 
0.241802E 
0. 6522 8E 


03 
03 
03 
03 


02 
02 
02 
02 


03 


134 


0.61 360€E 
0. 70610E 
0.79860E 
0.891165 


0.94806E 
0.60306E 
0,60306E 
0.85306E 


0.58656E 
0.67 O3SE 
02 f4886E 
0. 83528E 


0.90909E 
0.S53547E 
0 -45496E 
0.69967E 


0.59077E 
0.,65861E 
0.72427E 
0.f78108E 


..9278TE 

« 50223E 
0.43394E 
0.63720E 


03 
03 
03 
03 


02 
02 
02 
02 


03 
03 
03 
03 


02 
02 
02 
02 


03 
OS 
03 
03 


02 
02 
02 
02 
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oe Gen. 


a - or 


3 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0~2~54481E 03 
0.c0031E 03 
0.6G6781E 03 
O.273058E 03 
0.81607E 03 

VAPOR FLOW RATE 
O0223000F 02 0.94806E 
O.-91831E 02 0,89988E 
Oe-S5S1I724E 02 0.49854E 
0246192E 02 0.73103E 
0.67026E 02 


0 .56594E 
0.60842E 
0.68159E 
0.73849E 


03 
O03 
C3 
03 


02 
02 
02 
02 


O0.S7540E 
O.63176E 
Of7123253E 
O.74795E 


0.95160E 
O.53071E 
O0244503E 
0.74026E 


o SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0254474E 03 
0.60010E 03 
0.66798E O03 
O24 73050E 03 
6.81523E oO3 

VAPOR FLOW RATE 
O0e-23000E 02 0.94806E 
O0.89854E 02 0.838107E 
O.S1435EF 02 0.49707E 
Gea4ASISTE 02 O.WTPSTCE 
0.66316E 02 


0.S56S58SE 
0.60823E 
0.68176E 
OTST TSIE 


c3 
03 
03 
03 


C2 
C2 
C2 
02 


0.5 7533E 
0.63176E 
O27 13S7GE 
0.74E630E 


0.-95049E 
0.S201SE 
0244402E 
0.71929E 


0.46107E-03 


03 
03 
03 
03 


02 
02 
02 
02 


0.58282E 
0.6465 3E 
0.71 767E 
0.76050E 


0.94562E 
0.5311 0E 
0.43185E 
0.73907E 


0.14002E-04 


03 
03 
03 
03 


02 
02 
02 
02 


0. 58271 
0.64665E 
0.71 786E 
0. 75826E 


0.93999E 
0. SZ eISE 
0.2.42939E 
0. 7209S5E 


03 


03 
03 
03 
03 


G2 
02 
02 
02 
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~S91ISTE 
0.6S761£E 
0.72305E 
C.77964E 


0.93389E 
0.52609E 
0.44442F 
0.72266E 


O.S9141E 
0.65779E 
0.72322E 
0.77753E 


0.91960E 
0.52050€E 
0.44003E 
G8. 7DIDSE 


03 
03 
03 
03 


02 
02 
02 
02 


03 
03 
03 
03 


02 
02 
02 
02 
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ESE COMPUTED RESULTS ***%x% 
NUMBER OF ITERATIONS 4 
COLUMN PRESSURE (PSIA) 0.26470E 03 
REFLUX RATIC 3.212200 
PRESSUR » TEMPERATUR +» ENTHALFY AND LIQUID FRACTION OF FEED 
O02 265E 03 PSIA G2127E O32 F OF204E 06 BST ss. 0.0 
OUS2656 03 PSIA° G.2T1LE O39 F G.483E 06,.857Ts.0- 1,000 


xe STAGE VARIABLES *** 


PLATE NO. TEMPERATURE LIG. F. R. VAPS OF. R. HEAT DUTY 
7 LB MOL/YTIME LB MOL/TIME BTU/TIME 
1 0.85135E 02 0.71606F 02 0.23000F 02 -0.43341E 06 
2 0.1062SE 03 0.72049F 02 0.94806F 02 0.0 
3 0.11573E 03 0.70999F 02 0.95049F 02 0.0 
4 0.12311E 03 0,53960E 02 0.93999F 02 0.0 
5 O.13181E 03 0.51854E 02 0.91960F 02 0.0 
6 0.140S0EF 03 0.50107E 02 0.89854E 02 0.0 
7 0.14863E 03 0.4€515E 02 0.88107E 02 0.0 
8 O.172Z16E 03 0.48753E 02° .0.52015E 02 0.0 
9 O0.18705E 03 0.48550E 02 0.52253£ 02 0.0 
10 0.19819E C3 0.47935E 02 0.52050E 02 0.0 
11 0,20838E 03 0.46207E 02 0.51435F 02 0.0 
12 0,22216E C3 0.40902E 02 .0.4S707E 02 0.0 
13 0.25410E 03 0.10494E 03 0.44402E 02 0.0 
14 0.25826E 03 0.10600F 03 0.42939EF 02 0.0 
15 0.26362E 03 O0.1C7T1IS5SE 03 0.44003F 02 0.0 
16 0.27090E 03 0.10827E 03 O0.45151E 02 0.0 
v7 0.27813E 03 0.10893E 03 0.71270E 02 0.0 
18 0.28670E 03 0.1CS10E 03 0.71929F 02 0.0 
19 0.29866E 03 0.10811E 03 0.72095F 02 0.0 
20 0.31794E 03 0.10332E 03 0.71113E 02 0.0 
0.68 
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mh 
~~ 


21 0.35S563E 03 0.37000F 02 0.,66316E 02 
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SeintheQ®-  oMeeeel 
SIO THE LO ie oivee 


9.0 85 350666 .6 
0.0 9 #PE054,0 
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0.6 0 FOTSIT.9 
©.0 $0 S@keOty .0 
5,9 $0 S@vost.od 
O.0 GU ZEILIT.O 
‘SO 559£804,.0 G0 G6fE59.0 
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TOTAL 


Oe.S7OO00E 02 


xe MATERIAL BALANCES *x%% 
FEEDS 
PLATE NO. a PLATE NOs 13 
LB MOL/STIME MOL FRC, 1B MOL/STIME MOL FRC. 
METHANE 0.20000E O01 0.579 72-1 0.0 0,0 
ET HANE O-10000EF 02 C.28S986E 00 0.0 0.0 
PROPYLEN 0.-60000EFE O01 O.17391E 00 O.10000E O1 0.2526 75-01 
PROPANE O,12000EF C2 C.34783E CO OG. 7O00CCE O1 Dw.10687E GO 
IBUTANE O21O0000E O01 0. 28986E-01 0,.400C00E O01 0.6106CE-01 
NBUTANE O.30000E O1 0.869S57E-O1 O.17GO0E 02 0.25954E 00 
NPENT ANE O0.S0000E 090 0.144¢63F-01 0.215200EFE O02 O0.23206E 00 
HEXANE 0. 0.9 0.90000E O01 0.21374CE O00 
HEPTANE 0.0 0,0 Os.45000E 01 0,638702E-01 
OCT ANE 0.9 C20 0O.430C0E Cl 0.6564CEFE-01 
DECANE 0.0 020 0.35000E O1 0.2.53435E-01 
TOTAL 0.234500FE O02 0.65500E 02 
PRODUCTS 
VAPOR OISTILLATE IQGIO BiIStALeaArte 
LB MOL/STIME MOL FRC, iB MOL/TIME MOL FRC. 
METHANE 0.19599E O1 ¢. 652 13-01 0.0 0.210834E-01 
ETHANE 0.90252E 01 0.39240EFE 00 0,0 O.2i 275 9 
PROPYLEN O.38686E Ci C.1G6820E 09 0.0 0.22483E 00 
PROPANE 0.8013SE G1 0.34841E 00 0.0 0.53214E 00 
IBUTANE 0.88766EF- 01 0.38594E-02 0,0 0.12046F-O01 
NBUTANE 0.44060E-01 0. 1915S6E—-02 909.0 9.79 7286-02 
NPENTANE 0.94465E-05 0.41072E-06 0.0 0.4 2356E—-05 
HEXANE 0.30262E-10 G., PSRS fe— D2 0.0 0.34176E-10 
HEPTANE 0.20524E-14 6.89235E-16 0.0 0.52029E—-14 
OCTANE 0.13464E-18 0.58541E—-20 0.0 0,82294£-18 
DECANE 0.412879E-26 6. 1S208E-—2e 0.0 G,11835E-24 
TOTAL O.23000E 02 0.0 
BOTTCM FRODLUCT 
LB MOL/STIME MOL FRC. 
ME THANE 0.20328E-13 0254940E-15 
ETHANE 0.213228E-07 0. 357S50E-09 
PROPYLEN Ge27 IEZE“Ge 0.73466E-04 
PROPANE 0.41201&-01 0.111356-02 
IBUTANE 0, 70362E 060 6.1901 7E-01 
NPUTANE 0.52278E O01 0.14129E 00 
NPENTANE O.11285E 2 C.304S9E 00 
HEXANE O.79107E O1 O.2136€0E 00 
HEPTANE 0.42148E O1 O2.1L13S9iIE OO 
OCTANE 0.41540E O1 G.1Ll227E Go 
DECANE 0,.34605E 01 0.935 28E~-01 
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R000TE.O JATOT 


LIQUID SIDE 


ME THANE 
ETHANE 
PROPYLEN 
PROPANE 
IBUTANE 
NBUTANE 
NPENT ANE 
HEXANE 
HEPTANE 
OCTANE 
DECANE 
TOTAL 


STREAMS 
LIQUID ( 4) 
LB MOL/TIME MOL FRC. 


0.3991 7E-01 
O.97464E 00 
0.28738E 01 
0.87361E O01 
0.10082E O1 
0.213589E O1 
04.6465 0E-02 
0.81088E-06 
0.11570E-08 
0.20882E-11 
0.21706E-16 
O.15000E G2 


VAPOR SIDE STREAMS 


ME THANE 
ETHANE 
PROPYLEN 
PROPANE 
IBUTANE 
NBUTANE 
NPENTANE 
HEXANE 
HEPTANE 
OCTANE 
DECANE 
TOTAL 


VAPOR 
1B MOL/STIME 
0.29534E-08 
0.58430E-04 
0.25472E OC 
60.2207SE 01 
0.31985E 01 
0.135652 G2 
0.44093E O1 
O.10906EF Ci 
O.2S8S5T2Ee 00 
0.14641E O00 
0.5961 36-01 
O.25000E 02 


0.26612E-02 
C.64S 76E-01 
0.19159E 00 
0.58241E 00 
0567 22ie—08 
0.90590E-01 
C, 56433E-—03 
0.54059t- 07 
Ost TL SBOE—1TO 
0,139216—ic 
02-14470E-17 


(16) 

MOL FRC. 
0.11813E-09 
C.23372E-05 
C.101€SE-O1 
0.88299E-01 
C.12794E 00 
0.53470E 00 
0O.17637E 00 
C.43€25E-01 
0,114296-01 
0.58563E-02 
C.15845F-02 
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APPENDIX 3-8 TEST PROBLEM NUMBER - 8 


eee PROBLEM STATEMENT ***** 


TYPE OF COLUMN 3 

NO. OF PLATES 25 

NO. OF COMPONENTS 18 
PRESSURE AT TOP PLATE (PSIA) 


PRESSURE DROP PER PLATE (PSIA) 


NO. OF FEEDS 1 
ENTERING PLATE OF FEED 1 9 
FEED QUANTITY 0.58735E 03 
PRESSURE OF FEED 350.00 


TEMPERATURE OF FEED (F) 
REFLUX RATIO 290 


LIQUID RATIO OF TOP PRODUCT 


132 


359.00 


250.00 


0.5000 


SPECIFIED TOP PRODUCT 0.420408 03 


*** RATE GF CONVERGENCE WITH ITERATION NOW *** 


re) SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


02 57960E 
0.63897E 
0. 69835E 
O02. 7S772E 
O.81710E 


03 
03 
03 
03 
03 


O.5SG147E 
0.65085E 
0.71022E 
0. /76960E 
0.82897E 


VAPOR FLOW RATE 


O0221020E 
0.2.12192E 
0. 80035E 
0.80035E 
0-2 80035E 


03 
04 
03 
03 
03 


0.212192E 
9.220925 
0.80035E 
0.80035E 
0. 80035E 


C3 
O3 
C3 
03 
Ue 


04 
04 
03 
03 
03 


0.60335E 
0.66272E 
Oe-f2210E 
Ow.7EL4TE 
0.84085E 


0.212192E 
0.12192E 
0. &C03S5E 
0.80G35E 
0.80035E 


1 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


0.58149E 
0.69342E 
0. 70956E 
02.76189E 
Os8lLOSCE 


03 
03 
03 
03 
03 


0.61873 
0.69626E 
0.72008E 
O.77274E 
0.81926E 


VAPOR FLOW RATE 


Oe 21020E 
O02 17081E 
0. 80148E 
0.283934E 
0.88181E 


03 
04 
os 
03 
03 


O.12192E 
O.17818E 
0,81839E 
0.84668E 
0.88023E 


03 
03 
03 


0,63693€ 
02.66S972E 
O.73019E 
0s 183232 
0.82945E 


0.11€56€ 
0.14128E 
0.82601€ 
0.-85610E 
OC. SFT SE 


2 SUM OF SGUARES OF RESIDUALS 
TEMPERATURE ; 


0.S58420E 
0268718E 
02 73970E 
0.77939E 
0.281242E 


03 
03 


0.65030E 
0.69215E 
0.74725E 
0.78823E 
O.81701E 


VAPOR FLOW RATE 


0.21020E 
OwlV1IOTE 
02. 81 340E 
0.87337E 
0. 96189E 


03 
04 
03 
03 
03 


0. I12192E 
0217214E 
0. f4464E 
0.88751E 
0.96588E 


03 
03 
03 
03 
03 


C4 
04 
03 
03 
03 


0,.67C43E 
0.69386E 
0275430E 
0.7S9616E 
0-82406E 


0.12€C9E 
0.1 33562E 
0.854G91E 
0.90€73E 
0.95135E 


O.21987E 01 


03 
03 
03 
03 
03 


04 
04 
03 
03 
03 


0. 61522E 
0.67460E 
O.73397E 
0.79335E 
0. 85272E 


0.212192E 
0.12192E 
0. 8003S5E 
0.80035E 
0,80035E 


0.64382E 00 


O03 
03 
03 
03 
03 


04 
04 
03 
03 
03 


Ge 67324E 
0,68164E 
0.74045EF 
0. 79302E 
0.84351E 


0, 13652E 
0.12226E 
0. 823 029E 
0. 86677E 
0.84735E 


0.371 30E-01 


03 
03 
03 
Os 
03 


04 
04 
O03 
03 
03 


0.67614E 


0.706426 


0.7619 2E 
0.80278E 
0.8375 0E 


0, 15398E 
0,12205E 
0.85967E 
0.92 793E 
0.90839E 


03 
03 
03 
03 
03 


04 
04 
03 
03 
03 


03 
O02 
03 
03 
03 


O4 
O04 
03 
03 
03 


03 
C3 
03 
03 
03 


04 
04 
03 
03 
03 
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0.82710E 
0.68647E 
0.74585E 
0. S0S22E 
0.8&6460E 


0.212192E 
0,80035E 
0.,.80035E 
0.80035E 
0.8003S5E 


0.68710E 
0,€9726E 
O0.75104E 
0.80205E 
0.86390E 


0.15914E 
0.774S53E 
0.83420E 
0.87592E 
0.80939E 


0.68220E 
Oe 72838E 
0.77040E 
0.-80801E 
0.86531€E 


0.16645E 
0. 74187E 
0,80448E 
0.947S96E 
O.81346E 


03 
03 
03 
O3 
O03 


O4 
03 
03 
O3 
03 


03 
O03 
03 
03 
03 
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03 
03 
03 
Os 


03 
03 
03 
03 
03 
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03 
03 
03 
03 
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3 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0.58393E 03 
O268841E 03 
O275557E 03 
O.80271E 03 0.80780E 
0.81719EFE 03 0.82057E 

VAPOR FLOW RATE 
0e21020F 03 0,12192E 
O©e-12409F 04 0.12251E 
02.82299EFE 03 0.84948E 
0.88770E O03 0. 89655E 
O.90381EF 03 0.89818E 


O.,64761£F 
0. 69516E 
0.76648E 


04 
C4 
03 
C3 
C3 


0 .6€€6€40E 
0.70088E 
OL 77TETSE 
0.81106E 
0.82680E 


0.213368E 
0.12697E 
0.8€C079E 
O.90264E 
0,8&421E 


4 SUM OF SGUARES OF RESIDUALS 


TEMPERATURE 
0.58357E 03 
0.268803E 03 
0275544E 63 
0.2.79233E 03 O0.79737E 
G.8123S6E 03 6, 817TIE 

VAPOR FLOW RATE 
O02.21020E 03 0.12192E 
O.12951E 04 0,12697E 
O02 79665E 03 0.282007E 
0. 84079E 03 0.86044E 
GO. S9S253E 03 08. 9SITDGE 


0.64775E 
0.69380E 
0. 76467E 


03 
03 
03 
03 
03 


04 
04 
03 
03 
C3 


0.66754E 
O.70167E 
0.7 7257E 
0.280163E 
0.282445E 


0.212989E 
0.12419F 
OS825I1E 
0.88422E 
0.91 51V7E 


5 SUM OF SQUARES OF RESIDUALS 


TEMPERATURE 
0.583S7E 3 
0.68933E 03 
0.75726E 03 
0.79219E 03 0.79675E 
0.281167E O03 0.81655E 

VAPOR FLOW RATE 
O0.Z21020E 03 0. 22192E 
021302S5E 04 0.12821E 
O.82315E 03 0.86163E 
0293456F 03 0.94498E 
0.95204E 03 0.94342E 


0.64800E 
0.69539E 
0.766S52E 


C3 
03 
03 
03 
03 


04 
04 
03 
03 
03 


0.6E7E2E 
0270327E 
02.77422E 
0.-80073E 
0.82406E 


0.12009E 
0212492E 
0.&&453E 
0.9S157E 
0.92304E 


0.26399E-02 


03 
03 
03 
o3 
O03 


04 
O04 
03 
03 
03 


eee | 


03 
03 
03 
03 
03 


04 
04 
03 
03 
03 


0.67605E 
0.271528F 
0.78665E 
C.81320E 
0. 83921E 


0.12667E 
0. 12230E 
0. 86886E 
0.90666E 
0. 8506SE 


57 16E-03 


0.67670E 
0.71 592E 
0. 77S978E 
0.80527E 
0.83760E 


0. %t3213E 
0.11697E 
0. 82573£ 
0,90673E 
0. 87486E 


0215400E-04 


03 
03 
03 
03 
03 


O04 
O04 
O3 
03 
03 


0.67 735E 


Oef1721E 


0. 78098E 
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0.83 7355E 


0,23205E 
0.118436 
0. 90331E 
0.95465E 
0, 88164E 


14] 
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0O.f4166E 
0279556E 
0.81503E 
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0.&6263E 


0.213165E 
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0.95493E 
0.,80497E 
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NUMBER OF 


COLUMN PRESSURE 


REFLUX RATIO 


PRESSUR 


0. 350E O03 PSIA 


ITERATIONS 5 
(PSIA) 0.35000E 03 
1.¢°0000 
TEMPERATUR 
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ENTHALFY AND LIQUID FRACTIGN OF FEED 


**X* STAGE VARIABLES **% 


PLATE NCO. 


CONA MSW 


TEMPERATURE 


F 
0.12397E 
0. 18840E 
Os20822E 
O.227T 752 
O.22411F 
0.22973E 
0,23579E 
0.24367F 
0.25761E 
0.28440E 
0.29766E 
0.30693E 
0.31462E 
0.32138E 
0. 32735E 
0 733259E 
G.3SSTi1SE 
0, 341146 
0,34475E 
0.34826E 
0. 35207E 
0.35695E 
0, 36446E 
6377 T5E 
0.40303E 


Re 


LB MOL/STIME 


0.7S876E 
0.88054E 
C.SOO01LIE 
0.-89611E 
0,88212E 
0.86165E 
0.82879E 
0. 7E3S1E 
0.89972E 
0.9SO011E 
0. 10286E 
0.10515E 
0. 10703E 
0.10872E 
0.11015€ 
0.11119E 
0.2.11185E€ 
0.11216E 
0.11219E 
O0.11190E 
0. 11104E 
0. 10900E 
0.10486E 
0.971S62E 
0.16695E 


03 
03 
03 
03 
OFS 
03 
03 
03 
03 
03 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
O04 
04 
03 
03 


O.250E G3 F 10.5216 07 B.T..U. 


VAP. F. 


Re 


LB MOL/STIME 


0.221020E 
0,12192€E 
0. 13009E 
0.13205£ 
0, 13165E 
6. 150255 
GO. 12821 
0.212492E 
0.2118643E 
0,7327TE 
0. 82315E 
0.86163E 
0, 88453E 
0.90331E 
O.92021E 
0. S3456E 
0.94498F 
0.95157E 
0, 95465E 
0,.95493E 
0. 95204E 
0. 94342E 
0,92304E 
0. 88164E 
0.80497E 


03 
04 
04 
04 
O4 
64 
04 
04 
04 
03 
03 
03 
03 
03 
03 
03 
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%**e% MATERIAL BALANCES * 4% 


FEEDS 


cG2 
NITROGEN 
ME THANE 
ETHANE 
PROPANE 
IBUTANE 
NBUTANE 
IPENTANE 
NPENTANE 
2MPENTAN 
HEXANE 
HEPTANE 
OCTANE 
NGONANE 
DECANE 
UNDECANE 
DODECANE 
TRDECANE 
TOT AL 


PRODUCTS 


co2 
NITROGEN 
METHANE 
ET HANE 
PROPANE 
IBUTANE 
NBUTANE 
IPENTANE 
NPENTANE 
2MPENTAN 
HE XANE 
HEPTANE 
CCTANE 
NONANE 
DECANE 
UNDECANE 
DODECANE 
TRDECANE 
TOTAL 


PLATE NO. 


LB MOL/STIME 
0O2.41°900EF C0 
0.2L PCGE* OO 
O.74737E G2 
0.78890E 02 
0,13006E 03 
O0.44244E O02 
0.90173E C2 
0241344E O02 
0.43606E C2 
0215143E O2 
0.14955E 02 
O.25817E O2 
0.15875E O62 
0.87840E 061 
0.23920E O1 
0.44300E 00 
0.17400E 00 
0.88000E-01 
0.587355 03 


9 

MOL FRC. 

OQ. 7 S372—04 
0, 3S5924E-03 
0.12724E 00 
C.13431EF 00 
Q.22143EF 006 
QO. f53298E-—01 
0. 1S352E "00 
0.7039 1£-01 
C274242E-01 
0.257826~01 
0.25462E-01 
C.43955E-01 
0.27 028E-01 
0.14955E-61 
C.40725E- 02 
0. 7S423E-03 
C.29€24E-03 
0.14983E- 03 


VAPOR DISTILLATE 


LB MOL/TIME 
0.33030E 00 
0.19721E 00 
02.65482E C2 
0.52327E 02 
0.57040E 02 
0.12834E 02 
0021539E 02 
0.35464E 00 
0.94546F-01 
0.29464E-03 
0.79251F-04 
0.69443E-06 
0.15557£-08 
0.45104E-11 
0 9 3487E-14 
0.16997E-16 
0072871F-19 
0.44380F-21 
0.21020EF 03 


MOL FRC. 
0. 15713&-02 
0.93818E-03 
0.31152E 00 


C.24ES4E 00. 


0.27126£ 00 
C.61057E-01 
C210247E 00 
0. 16872E-02 
0.4497SE-03 
0.14017E-05 
0.37703E-06 
C.33027E-08 
0e74009F-11 
0.21458E-13 
0244475F- 16 
0,80863F-19 
C. 34667E-21 
0.211136-23 
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PLATE NC. 


LB MOL/TIME 


MOL FRC, 


LIQUID DISTILLATE 


LB MOL/STIME 
C. BB 7S0E-01 
0.13987E- 01 
0.92990EFE O1 
0.2656SE 02 


-0.72996E 02 


0,31302E 02 
0.67197E 02 
0.20552F O01 
0.67307E 00 
0.39712E-02 
0.12820E-02 
0. 23475E- 04 
6211703E-06 
C. 70652E-09 
0. 2886SE-11 
0.98804E-14 
0. 7T8898F-16 
0.87851E-18 
0.21020E 03 


MOL FRC. 
02.42241E-03 
0.6654 3E-04 
0.4423SE-01 
0.1264C0E 00 
0.34727E O00 
O0.148S1E 00 
0e-31968EF 00 
0.97772E-02 
0.0202 16-02 
0.18892E-04 
0,.60990E-05 
0.11168E-06 
0.5567 4E-09 
Gece lee—it 
G, 137 32E-14 
0.4700 56-16 
0.3 7S3S5E-i19 
0.41794E-290 
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&§-3Veeto. 0 
Al-a1T2041.0 
&2-41T8ST60 
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cGe2 
NITROGEN 
METHANE 
ETHANE 
PROPANE 
IBUTANE 
NBUTANE 
IPENTANE 
NPENTANE 
Z2MPENTAN 
HEXANE 
HEPTANE 
OCTANE 
NONANE 
DECANE 
UNDECANE 
DOCDECANE 
TRDBDECANE 
TOT AL 


BOTTCM 
iB MOL/TIME 
0. 160563E—-12 
0.38112E-16 
0.42 723E-12 
0.11606E-08 
0,24584E-03 
0.98097E-01 
O.2i417V2E 01 
0. 38933E 02 
0.42838E 02 
0.1 S51S9E O2 
O.14954E O2 
9.258175 62 
0.15875E Ce 
0O.87840E Oi 
0.23920E O01 
0.44300E O00 
O217400E O00 
0.88000E-01 
0.216695E 03 


PRODUCT 


MOL FRC, 
0.99 8075-15 
0.22828E-138 
Ce 25S5S0E-14 
0.69518E-—11 
G.14725E—-C05 
0. S87S8E— 03 
0.84887E-O02 
8.23320& 0 
0.25659E 06 
0.90677E-01 
0. 89S€68E- 01 
0.21546€4E 00 
G. 95 087E-O01 
C.526€14E-01 
0,14327E-01 
5. 265SSE-02 
0.210422E-02 
0,52719E-03 
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APPENDIX 3-9 TEST PROBLEM NUMBER - 9 


*x*kxe PROBLEM STATEMENT **% xx 


TYPE OF COLUMN 3 
NO. OF PLATES 70 
NO. OF COMPONENTS 4 
PRESSURE AT TOP PLATE (PSIA) 84.00 
PRESSURE DROP PER PLATE (PSTIA) 0.0 
NO. OF FEEDS 1 
ENTERING PLATE OF FEED 1 35 
FEED QUANTITY 0.18195E 04 
PRESSURE OF FEED 84.00 
LIQUID RATIO OF FEED 1.0000 
REFLOX RATIO 11.300 
LIQUID RATIO OF TOP PRODUCT 0.0 
SPECIFIED TOP PRODUCT 0.67000EF 03 


25.31 1.039) RC=0 


0 SUM OF SQUARES OF RESIDUALS 
TEMPERATURE 


O0.57160E 
O.». STAS4E 
0. S7508E 
0.5768 1E 
0.257855E 
0.58029E 
0.58203E 
0258377E 
0.-58551E 
0. 58724E 
0. 58898E 
0259072E 
Oe 59246E 
0259420E 


03 
03 
03 
03 
03 
03 
OS 
03 
03 
03 
03 
03 
03 
03 


O.571S5E 
0.2.57369E 
0.57542E 
0.57716E 
0.57890E 
0.S8064E 


S«SIZEVE 
0.59455E 


VAPCR FLOW RATE 


0267000E 
0. 82410E 
O02 S2410E 
0.82410E 
0282410E 
0.282410E 
0.8241 0E 
0282410E 
0282410F 
0. 82410E 
0. B2410E 
0.82410E 
0.82410E 
0.82410E 


03 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 


0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410EF 
0.82410E 
0.82410E 
0.,82410E 


03 
03 
03 
03 
03 
03 
03 
03 
C3 
03 
03 
0S 
03 
C3 


04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 


0.S57230E 
O.57403E 
O.S7577E 
Go .S77St 
0.5792S5E 
O.5€CS9E 
0.S58272E 
Oe. S8446E 
0.S€620E 
0 .S8794E 
0.58S68E 
0259142E 
O.S9SPSE 
0.5S489E 


0.82410E 
0.82410E 
0e.62410E 
0.82410E£ 
0.82410E 
0.82410E 
0.824106 
0-82410E 
0282410E 
0.82410EF 
0.82410E 
0.82410E 
0.82410E 
0.-82410E 


*%** RATE OF CONVERGENCE WITH ITERATION NO. 


KE He 


0.92701E-02 


03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
O04 
04 
04 
04 
04 
04 
O04 
04 
04 


0.57264E 
0,57438E 
0.57612E 
0. S7786E 
0.57 960E 
0.58133E 
6,.58307E 
0.58481E 
0. S8655E 
0.2.58829E 
0.59003E 
0. S59 176E 
0.59350£ 
0.59524E 


0.82410E 
0.82410E 
0.82410E 
0.82410F 
0.82410E 
0.82410E 
0.8241 0E 
0.82410E 
0.62410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 


03 
03 
03 
03 
03 
03 
03 
03 
03 
os 
03 
03 
03 
03 


04 
C4 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
O04 
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0.-57299E 
0.57473E 
0.57647E 
0O.57821E 
0.57994E 
0.58168E 
0.258342E 
0.S58516E£ 
0.58690E 
0.S8864E 
0.>590S7E 
0.5921 
0. SISSS5E 
0. 59560E 


0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 
0.282410E 
0.82410E 
0.32410E 
0.82410E 
0.82410E 
0.82410E 
0.82410E 


03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 


04 
04 
O4 
04 
04 
04 
04 
04 
04 
04 
04 
O4 
04 
04 


Wecyve.t) 
3etec2.0 
3v#atvTe. 0 
afeeyvye.o0 
sa0CT2.6 
346 182.0 
J546E42,.0 
$o1256¢.,.0 
#0° 4582 .0 
364 Oe 2. Oo 
ITEC? ,.O 
sj 7re¢e.0 
,2n€ o2.0 
200702 .0 


207 4%6.0 
90 1455.0 
791058,.0 


19 '658,.0- 
201 450,0 | 


1OLece,¢O 
Of218.0 
40169A.0 
101 ace.O0 
90 §@98 .0 
9Ofeso,6 
30148 ,@ 
3O/e94.G 
20 £488.0 


oe .0OM WT 


SO-210T¢SC.0 


SeoctT2.0 


SREP Te.O 


aState.o 
S6aX $2.0 
S040 T2.0 
Ste (62 .O 
S* Oot 62,0 
(ene .O 
$2294 .0 
ave 82.0 
2.0004 .0 
253 1069 .0 
sOCRGF 0 
wie! 0 


901466,.6 
40 1658.0 


orese.o0 


30!:46386,.0 
3 16S4.0 
3614¢@¢S6,.9 
90 16¢64,0 
402+ Ge .0 
30 1456.0 
wies4.0 
40 1*24,0 
701658,0 
goresa.o 
79 14¢8.0 


904 30 14464.0 CO oo 

2s tnanagto 4016886 -0 £9" 0ggtane 

30 PO 701499.0 #0 BOT8SB.O 
90 14S8.0 


oo S0T0S4.0 80 a01740%,.0 40 

of 30TO6@.0 60 SOrece. 0 #6 FOTSS8.0 
eo 3018S 20 701856,.0 80 3O14#88.0 
“0 30(0S4.6 60 301444.0 “0 2076580 
AO 301466 50166°.5 © aGreS6.0 
60 308 SO6ela.e a0 SO1aee co 
#0 40fe OO 301458.0 #0 S01 888 
60 SOb0EGS0 20 GOTSSS,O 60 seounee 
*0 408 8 301e5n.0 66 901894.0 
oO MOLE #0 301456.0 © 961488,0 
*0O SOOTHES OG 40 FO14SH.0 80 3010S8,.0 


1 SUM OF SQUARES GF RESIDUALS 


0.2.57282E 
0.2 57368E 
0.257432E 
0.5757 2E 
04 5782S5E 
Owls ee 
0.. 8850 3E 

a>O?r23e 
0. 5S8904E 
0.59129E 
0. 59353E 
0.59514E 
0.2 59586E 
0. S9600E 


TEMPERATURE 


03 
03 
03 
03 
O03 
03 
03 
03 
03 
03 
03 
03 
03 
03 


Os STS23E 
O, STSTISE 
0.57452E 
0.57612E 
0.257888E 
0. 58250E 
0.58SS57E 
0. 58756E 
0.-58947E 
O~59 17 Te 
0.59393E 
0.59535€ 
0.59592E 
0.59602E 


VAPOR FLOW RATE 


0.267000E 
0. 82392E 
0.82286E 
0. 82025E 
0. 81545E 
0-2 SO894E 
0. 80205E 
Ow 79696E 
0. 79624E 
0.79451E 
O.79191E 
Oe 78847E 
0. 78460E 
0.78108E 


03 
O04 
04 
04 
04 
O4 
04 
04 
O04 
04 
04 
04 
O04 
04 


0.82410E 
0.82405E 
0.82249E 
0281944E 
0.81425E 
0.80727E 
0.80073E 
0.796951E 
O.TIS96E 
0. 79406E 
0.79126E 
0. TST 70E 
0. 78382E 
0..78040E 


O03 
03 
03 
C3 
03 
C3 
03 
03 
03 
03 
03 
C3 
03 
03 


04 
04 
04 
o4 
C4 
04 
C4 
04 
04 
C4 
04 
04 
04 
04 


0 .Sfa33e 
0.S57391E 
0,57476E 
0. S785S9E 
O.57S61E 
0,.S8S317E 
0,.5€€04E 
0.58798E 
0.58990E 
0.59223E 
0259428E 
0.5 SS352E 
0.595966 
0.59606E 


0.82438E 
0.82354E 
0282203E 
0.8 P8525 
0.81283E 
0.80619E 
0,79949E 
0. /7S681E 
0.79566E 
0.79361EF 
0.79064E 
0.278694E 
O45 7TEZIGE 
OSTTSaCTE 


0.65545c- 04 


03 
03 
O03 
03 
03 
03 
03 
03 
O3 
03 
03 
03 
03 
03 


04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 
04 


O.57341E 
0. 573985 
6 PRS fa he re 
0. 57 708E 
0.58027E 
O-_ DESO 
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RAK MATERIAL BALANCES *%* 


recosS 
PLATE NG, 3S PLATE NC. 

LB MOL/STIME MOL FRC. LEB MOL/TIME MOL FRC. 
PROPANE 0.10000E O02 C.5496€0E- 02 
IBUTANE 0.65200E 03 0.258345 OC 
NBUTANE 02114006 04 0.62655E 00 
IPENTANE 0.17500E C2 0. 961 E0E— 02 
TOTAL 0,18195E 04 
PRODUCTS 

VAPOR DISTILLATE LIQUID DISTILLATE 
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